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O ' Abstract 

£T) | The oscillating gravitational field of an oscillaton of finite mass M causes 

it to lose energy by emitting classical scalar field waves, but at a rate that is 
non-perturbatively tiny for small fi = GMm/hc, where m is the scalar field 
mass: dM/dt ps -3797437.776333015(c 3 /G)^ 2 e- 39 - 433795197160163 ^[l + 
j> \ 0(n)]. Oscillatons also decay by the quantum process of the annihilation of 

•rH , 

^ , scalarons into gravitons, which is only perturbatively small in /i, giving by 

^ itself dM/dt « -0.008 513 223 934 732 692(m 2 c 2 /ft),u 5 [l + 0(p 2 )]. Thus the 

quantum decay is faster than the classical one for /j, ~ 39.4338/ [In (hc/Gm 2 ) + 
7 In (1/h) + 19.9160]. The time for an oscillaton to decay away completely into 
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free scalarons and gravitons is idecay ~ 27i 6 c 3 /G 5 m 11 ~ 10 324 yr(l meV/mc 
Oscillatons of more than one real scalar field of the same mass generically 
asymptotically approach a static- geometry U(l) boson star configuration with 
H = Ho, at the rate d(GM/c 3 )/dt « \(Cj^)e~ a ^ + Q(m/m Pl ) 2 h 3 ](h 2 ~ Mo), 
with Ho depending on the magnitudes and relative phases of the oscillating 
fields, and with the same constants C, a, and Q given numerically above for 
the single-field case that is equivalent to Ho = 0. 
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1 Introduction 



Seidel and Suen QjJ |2] have found numerically that there exist non-singular oscillat- 
ing self-gravitating solitonic configurations of a real scalar field, which they called 
oscillating soliton stars. These have also been studied by several other authors 
[HI |U 13 El M an d are now generally called oscillatons. In the simplest case, 
which is what I shall consider here, they arise from the Einstein-Klein-Gordon (EKG) 
equations for gravity plus one or more minimally coupled massive real scalar fields. 

The previous numerical evidence suggested that although these oscillatons are 
oscillating, they appeared to be periodic and stable [IJI21II100II10III], so that 
classically, at least, an isolated oscillaton might be expected to last forever. How- 
ever, here I shall show that oscillatons of finite mass actually decay classically as 
the oscillating gravitational field leads to the emission of scalar waves. The decay 
rate is calculated for the case of low-mass classical oscillatons and is found to be 
nonperturbatively tiny (nonanalytic in the oscillaton mass at zero mass), given by 
Eq. (j!35|) below, which may be why it has not been clearly seen numerically. 

Seidel and Suen [I] did recognize that their numerical results were consistent with 
quasiperiodic oscillations analogous to the orbit of two black holes that spiral inward 
while emitting gravitational waves. My results are also similar to this analogy, 
with the oscillatons classically emitting scalar waves instead of gravitational waves, 
except that here the classical decay rate goes to zero faster than any power law of 
the appropriate small parameter (here //, which is the mass M of the oscillaton, 
multiplied by the scalar field mass m, and divided by the square of the Planck mass 
mpi) as fi is taken to zero. 

In addition, oscillatons decay quantum mechanically by the annihilation of scalarons 
into gravitons, at a rate that is also calculated here for low-mass oscillatons, given 
by Eq. (j!52|) below. Although this mass-loss rate is also small, it is perturbative 
and goes as the fifth power of the oscillaton mass, so for sufficiently small oscillaton 
mass, this quantum decay dominates over the classical mass-loss rate. The time for 
an oscillaton with an initially large number of scalarons to decay away completely 
into free scalarons and gravitons then goes as the inverse 11th power of the scalaron 
mass and hence is very large if the scalaron mass is much less than the Planck mass. 

Although my numerical results for the classical and quantum decay rates are 
for a single real scalar field in a spherical configuration (and for any number of 
such fields of the same mass oscillating in the same mode except for possible phase 
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shifts [7j, which do give a nontrivial effect), I shall start by giving the formalism 
for the classical decay rates for an arbitrary nearly-Newtonian configuration of an 
arbitrary number of massive scalar fields, and then do a detailed numerical analysis 
of the single-field nearly-Newtonian and nearly-periodic spherical case for both the 
classical and quantum decays. Then I shall return to a discussion of the classical 
and quantum decay rates for multi-field oscillatons. 

Throughout this paper I shall assume that the mass of each scalar field is much 
less than the Planck mass, which is a necessary (though not sufficient) requirement 
for doing a classical analysis and is also necessary for the validity of various equations 

1 shall use for the quantum decay of oscillatons. 

2 Notation and Units 

Consider the case in which there are n real scalar fields each with mass mi, 
minimally coupled to Einstein gravity, and with no other self-interactions. The index 
/ labels the different mass values, and the index J labels the different scalar fields 
that have the same mass. 

In the classical analysis, I shall often use units in which c = 1 (though sometimes 
for results I shall insert the appropriate power of c in order to be able to evaluate 
quantities in conventional units), but I shall not set h or G equal to unity. However, 
to avoid having fr's in most of my equations, I shall let the masses mi have units of 
inverse time, which is indeed what they would have (at least if c = 1) in the classical 
Klein-Gordon equation that each scalar field obeys, 



I.e., a free zero-spatial-momentum real scalar solution in flat Minkowski space- 
time in orthonormal Minkowski coordinates, with the appropriate zero of time, 
would have the time dependence cos (m/t) . In terms of the conventional scalaron 
particle masses, which I shall hereafter denote with the starred subscript, m*i, one 
has mi = m^ic 2 /%. Perhaps it is more perspicuous to write this relation as 



so that in terms of the classical quantity mi (the natural frequency of the scalar 
field, in radians per second), the energy m*/c 2 of a one-particle quantum excitation 
of the scalar field is indeed h times the frequency of the excitation. 



(a-mj)$ij = 0. 



(1) 



m ^c 2 = hmi, 



(2) 
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That is, I am taking the view that it is the natural frequencies mj that are the 
classical parameters of the scalar fields, and that the masses m*j of the scalaron 
particles are quantum properties that will not show up in a classical analysis or 
in the classical decay of the oscillatons (though they will when one considers the 
quantum annihilation of scalarons into gravitons). 

Analogously, to avoid factors of Newton's gravitational constant G in most of my 
equations, when I consider the gravitational mass M of a scalar field configuration 
or oscillaton, it is convenient to include G in it (or actually G/c 3 if one uses units 
in which the speed of light, c, is not unity), so that my M has units of time and 
is thus half the gravitational (Schwarzschild) radius of the configuration divided by 
the speed of light. Therefore, if I let M* be the mass in conventional mass units 
(e.g., grams or kilograms), what I shall use is 

M = —5—- (3) 

With these conventions, I can avoid using h and G in most of my intermediate 
equations, even without using units in which those quantities are set equal to unity. 

For example, the simplest classical spherical oscillatons of a single real scalar field 
(with no nodes) are characterized (up to the overall scale, into which the natural 
frequency m enters) by the single dimensionless parameter 

GM*m* n / M* \ / m*c 2 \ , , 

fi = Mm = — ^ « 7.483 138 84 x 10 9 — - , 4 

he \M Q ) \ eV )' y> 

where M & pa 1.989 x 10 33 g is the mass of the sun in conventional units. (In temporal 
units, the mass of the sun is 4.925 490 95 x 10 -6 s, almost 5 microseconds, known to 
much higher accuracy than in conventional units, since the gravitational effect of the 
sun, proportional to GM Q , is known much more accurately than G is in conventional 
units.) By dimensional analysis, one can then easily see that if there is a classical 
decay of such an oscillaton, one must have (since my M has the dimension of time) 

£~/oo. 00 

a function purely of the only dimensionless parameter of the oscillaton, /1. (In Eq. 
(I135J) below I shall give this function for /x «C 1, finding that it is nonanalytic at 
/i = 0.) 

The price of this simplicity in the units is that one must get used to the mass of 
the oscillaton having the dimension of time, which is the inverse of the dimension of 



4 



frequency that is the classical 'mass' of the scalar field in its classical Klein-Gordon 
equation. 

We can get a further simplification by using an appropriate redefinition of the 
scalar fields Since the square of the time derivative of a scalar field has the 

dimension of energy density, the square of a scalar field has the dimension of mass 
divided by length, which is the same as the dimension of c 2 /G. Thus (considering 
also the 8tt in Einstein's equations) it is convenient to define the dimensionless scalar 
field values 

faj^^SnG/c^u. (6) 

Then by Einstein's equations, the Ricci tensor generated by the stress-energy 
tensor of the scalar fields is 

R »P = X][07J,a0/J,/3 + ^9apmi 2 (p 2 Ij]- ( 7 ) 

ij 1 

This will have the dimension of inverse time squared if the coordinates have the 
dimension of time and if the metric components g a p are dimensionless. 

Following the examples of [TQ|E|, it is also convenient to combine each rescaled 
dimensionless real scalar field (fiu and its time derivative <pu = dcpu/dt into a single 
dimensionless complex quantity, 

y IJ = l e im *\ci )I J + —hj)i (8) 
2 rrij 

so 

4>u = y u e- mit + y u e imit (9) 

and 

fa = -im^ue' 1 ™ 11 + zm/^ 7 je im/ *. (10) 

In terms of the complex tyjj and its complex conjugate ^?u, the time-time, time- 
space, and space-space components of the Ricci tensor are (using to denote the 
time coordinate t = x° and lower-case Latin letters to denote spatial coordinates x l ) 

Roo = EW?[( 2 + flbo)*/J*« " (1 " lgoo)^jje- 2tmit + ^jje- 2imit )]}, (11) 
i j 1 

Roi = Y,{ m h^iJ^iJ + \^)je- 2imit + ^ 2 ije 2imit )} 
u 1 

+ tmj^u^jj^ - *u% m - ^u^jj.e- 21 ^ + ^jj^jj.e 2 ^' 1 }}, (12) 
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ij / 
+ Vij^ij,! + *ij,i*uj + *u,i*ijjt + ^u,^iJ,je 2imit }. (13) 

3 Gauge or Coordinate Conditions 

In finding solutions to the Einstein-Klein-Gordon equations, one must make a choice 
of coordinates or gauge for the gravitational field. I shall restrict consideration to 
3+1 dimensional spacetime. There there are four coordinates to be chosen, giving 
the freedom of four free functions over spacetime for the gauge group of coordinate 
t r ansf or mat ions . 

I generally find it convenient to use these four degrees of freedom to set the 
time-space components of the metric to be zero, g 0i = for the three i, and to 
set g o to be independent, or nearly independent, of the time coordinate t. This 
then implies that the hypersurfaces of constant t are orthogonal to the worldlines of 
constant spatial coordinates x\ and that along each such worldline, the proper time 
is nearly proportional to the coordinate time t (with a space-dependent constant of 
proportionality). For example, if the metric is periodic in time, we can choose g 00 
to be precisely independent of the time coordinate t. However, this still leaves the 
freedom to make arbitrary spatial coordinate transformations that are independent 
of t. 

If we wish to pin down the spatial coordinates, we could, for example, choose 
them so that the time average of the spatial metric, (gij), over a time that is long 
with respect to the reciprocal of the smallest natural frequency difference \mi — mp\, 
is as nearly proportional to the 3x3 identity matrix as possible. More explicitly, 
if we take (g^) to be a spatially dependent 3x3 matrix, we could choose spatial 
coordinates so that the integral over all space of the square of the traceless part of 
this matrix is minimized. 

This would still not pin down the spatial origin or angular orientation of the 
resulting quasi- Cartesian spatial coordinate system, but we could choose the spatial 
origin to be that which gives the center of mass of the asymptotic form of the time- 
averaged metric. If it is necessary to fix the orientation, we could, for example, 
fix it so that the asymptotic quadrupole moment has its principle axes lying along 
the three coordinate axes in some order determined by, say, the ordering of the 
eigenvalues of the quadrupole moment. Of course, this specification is degenerate in 
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the spherically symmetric case, but then the angular orientation about the center 
of mass makes no difference. 

For our purposes below it is not necessary to be so precise, but I am just illus- 
trating how for a generic nearly periodic metric, it appears to be possible to fix all of 
the coordinates completely. Of course, there is some arbitrariness in the procedure 
chosen for this (e.g., whether to take the time average of the spatial metric or of 
its inverse or of some other matrix function of the spatial metric, and how to define 
the time average over a finite time if the non-oscillatory part of the metric is slowly 
varying). But once a sufficiently precise procedure is chosen, the coordinates are in 
principle rigidly given, and hence so are the metric components for a given spacetime 
in that coordinate system. That is, the procedure makes the resulting coordinates 
and metric components become procedure-dependent but gauge-invariant functions 
over the spacetime. 

As a result of all but the last parts of the procedure outlined above, one can 
write a generic periodic or approximately periodic metric in the form 

ds 2 = -e 2U{xk) dt 2 + e -m* k )+mx k ){[i + 2W(t,x k )]5 l] + a l3 (x k ) + h^t.x^dx'dx^ 

(14) 

where the time averages of the time-dependent quantities, that is the scalar W(t, x k ) 
and the traceless symmetric tensor hij(t,x k ), are all zero, and where the spatial 
coordinates are chosen to minimize the integral over all space of the trace of the 
square of the time- independent traceless symmetric tensor Uij{x k ). 

(If the metric is only approximately periodic, it may be able to be written exactly 
in the form above for only a limited amount of time. Alternatively, to be applicable 
for longer times, either the form above may be only approximate, or else one might 
need to give U, V, and some slow time dependence to deal with slow nonperiodic 
changes in the geometry.) 

Of course, there are many other similar forms for the metric from similar pro- 
cedures, such as having W(t,x k ), hij(t,x k ), and/or a i j(x k ) ) or suitable multiples of 
these quantities, as arguments of an exponential, so I am not claiming that there is 
a unique preferred form for a periodic or approximately periodic metric, but only 
that the form above, or its slight generalization to the case when U, V, and have 
some slow temporal variation, is sufficient for our purposes. 

Note that in the spherically symmetric case, to which I shall turn later, the time- 
averaged spatial metric is necessarily conformally flat by its spherical symmetry, so 
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(Tij(x k ) = 0. Both time-dependent quantities, W(t,x k ) and hij(t,x k ), are generically 
nonzero. The spatial dependence of any scalar quantity is a function of the one 
spatial function r 2 = 5ijX l x^ so U and V are functions purely of the 'radius' r, 
and W is a function of t and of r. With spherical symmetry, the traceless tensorial 
quantity has the form hijdx l dx :j = h(t, r){5ijdx l dx^ — 3dr 2 ) for some function h(t, r) 
of both time and radius whose time average is zero. 

4 Nearly-Newtonian Scalar Field Configurations 

In this paper I shall focus on self-gravitating configurations of one or more massive 
scalar fields in which the gravitational field is very weak (given to an adequate ap- 
proximation by the linearized Einstein equations), and the scalar fields have a very 
slow spatial dependence (so the dominant piece of ^fjj has a very slow spacetime 
dependence, though <pu does have a temporal oscillation of frequency roughly mj 
that is not considered slow, since slowness is taken to be relative to these frequen- 
cies). See |1U[ l§] for previous analyses in this limit, which have been a motivation 
for some of my choices of variables. 

In this limiting case, the metric functions U(x ) and W(t,x k ) are much smaller 
in magnitude than unity (but not negligible), and V(x k ), (Jij(x k ), and hij(t,x ) are 
negligibly small. Therefore, the metric (fPfj) takes the form 

ds 2 « -[1 + 2U{x k )]dt 2 + [1 - 2U(x k ) + 2W(t, x k )]5 ij dx i dx j . (15) 

If W(t,x k ) = 0, then the approximate metric (|15|1 would be truly Newtonian, 
but the temporal oscillations of the scalar fields give oscillating components of their 
stress-energy tensor and hence of the Ricci tensor components and of the 

metric (mainly at twice the frequencies of the fields themselves), so W(t,x k ) is 
nonzero even at the linearized gravity level. 

As we shall see below, since the scalar fields are assumed to have slow spatial 
dependences, a typical magnitude of W(t, x k ) (say its rms value at some spatial 
location where that is maximized, but of course not its time average, which is zero 
at all spatial locations by definition) is of the same order of magnitude as a typical 
magnitude of the square of U(x k ), which we have dropped in expanding the expo- 
nential. It is also of the same order of magnitude as a typical value of V(x k ), which 
we have also dropped. Therefore, it may be thought a bit of a cheat to include the 
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W(t, x k ) term in the metric above but not the U(x k ) 2 and V(x k ) terms, which are 
similar in magnitude. 

However, the point is that W(t, x k ) is the largest time-dependent term and is 
responsible for the dominant contribution to the classical decay of oscillatons and of 
other self-gravitating real scalar field configurations in the nearly-Newtonian limit. 
The U(x k ) 2 and V(x k ) terms that have been dropped are smaller than the U(x k ) 
time-independent terms that have been kept, and none of those terms directly con- 
tributes to the classical decay. Thus the philosophy is that the metric (fT5j) includes 
the dominant time-independent corrections to the flat Minkowski metric (the U(x k ) 
terms) and the dominant time- dependent corrections to the flat metric (the W(t, x k ) 
term, in the gauge in which g 0Q is independent of t by construction). 

Before going to periodic configurations (in the approximation of neglecting the 
scalar field emission), let us consider the slight generalization to the metric 

ds 2 « -(1 + 2U)dt 2 + (1 - 2U + 2W)5 ij dx i dx j (16) 

in which W has a time dependence at frequencies that are roughly twice that of the 
mi, but U is now allowed to have some time dependence that is even much slower 
than its slow spatial dependence. 

Now, instead of Eqs. ©-(HH, I shall take 

4>u « i[>ue- imit + ^je imit (17) 

without the restriction ifjjje~ irnit + x/j IJ e' imit = that is true for from Eqs. (|5)l- 
(110)1 . Instead, I shall assume that each 0/j is such that ipu can be chosen to give 
4>u approximately and also give 

\$u\ <mj|^jj| <t: mj\ipu\. (18) 

Then the Klein-Gordon equation in the metric (jlfi)) with \W\ \U\ 1, which 

is 

4> u w -mfyu - 2mjUfaj + c 2 V 2 <p u (19) 

when for the moment we ignore the W term, implies that each ipu approximately 
obeys the Schrodinger equation 

ic 2 

ipu « ^V 2 ^/j - inuUfaj, (20) 
2mi 

where V 2 is the flat-space Laplacian, V 2 ?/>/j = ^ipu^j. 
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Note that it is m*ic 2 U = HmiU that is the Newtonian potential energy of the 
particle of 'mass' mi, and not U itself (which is dimensionless) . Also, by having mi 
have units of frequency rather than conventional mass units, the explicit appearance 
of H is avoided in Eq. (|20p. This is what one would expect, since this Schrodinger 
equation came from the purely classical Klein-Gordon equation for the real scalar 
field, rather than from any quantum equation. Note that I have chosen ipu to be 
dimensionless rather than, say, having the spatial integral of its absolute square be 
unity (or perhaps some other positive integer), as one would normally normalize the 
wavefunction of a truly quantum Schrodinger equation. 

It is most straightforward to regard the approximate equivalence between the 
real second-order Klein-Gordon equation and the complex first-order Schrodinger 
equation (|2*U|) as a procedure that works when one starts with a solution of the 
Schrodinger equation (f2T)|) (with a weak gravitational potential, \U\ <C 1) for which 
|V 2 ?/'/j| <C mj\ipu\ (except near possible zeros of ipjj) and then uses Eq. (f^Uj) to 
construct from it an approximate solution of the Klein-Gordon equation. 

In the reverse direction it is a bit more subtle. If one uses Eq. (JHJ) to define \1/ u in 
terms of <pu and its time derivative, this itself will be close to the solution 
ipu of the Schrodinger equation (|2U|). However, since ^ u will generically have a 
small term roughly proportional to e 2tmit as well as its dominant term with a much 
slower time variation, the time derivative ^ u will pick up a relatively significant 
contribution from the term that is roughly proportional to e 2imit and so will be 
significantly different from ipu. Thus \1/ U defined by Eq. (JHJ) will not satisfy the 
Schrodinger equation (J2*njl . 

However, one can instead define 

i ■ 1 ... 

= Vu + - — *n = —^e mit (mj0 u + 2im I <t>u - faj) 
2mi Arrij 

~ y mi \<Pu+— iu + Ufaj-^V^j), (21) 

where for the last expression I have used the approximate form (|19|) of the Klein- 
Gordon equation in the Newtonian part of the metric to evaluate the second time 
derivative of <pu in terms of its value and its spatial Laplacian. This ipu then obeys 
the Schrodinger equation (|2U|) when U is small and slowly varying and when 0/j 
is oscillating at nearly its natural frequency mi and has a slow spatial variation in 
units of mj. 

With this definition of the complex ipu in terms of the real <pu and its derivatives, 
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Eq. (JT7j) is still a fairly good approximation for <pu in terms of ipu, but an inversion 
of Eq. (|21jl that is accurate to one higher order is 



TO; VTLj 

c 2 - r 2 



(1 - U)[(^j + ^ w )e^ + + ^ w )^. (22) 
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Nevertheless, although this formula for 0/j is a more accurate inversion of Eq. (|2~T|) 
than is Eq. 1)170. I don't know that it really gives a more accurate solution of 
the Klein-Gordon equation than Eq. ()T7)l does from a solution of the Schrodinger 
equation l)2U|) . 

After getting an approximate solution of the Klein-Gordon equation in the 
nearly-Newtonian metric (jTS)) with \W\ <C |Z7| <C 1 (and temporarily ignoring the 
small effect of the tiny but rapidly-time-varying W term, which shall be discussed 
later), we need to solve the Einstein equation for the effect of the stress-energy ten- 
sor of the scalar fields on the metric. To leading order, the resulting Ricci-tensor 
components from the stress-energy tensor and the Einstein equation are 

^oo « - imf^j} « EW/^w - l^je- 2 *™'* + $je- 2imit )}}, (23) 

R 0l « 0, (24) 

J2y « E ^iM?J « 5>fo[l^« + k^jc-*"'* + ^e 2 ^*)]. (25) 
u z ij z 

The corresponding Einstein tensor components are simpler, 

Goo«E 2m /^"> ( 26 ) 

G 0i « 0, (27) 

Cii « -E^WjC-^ + ^jC*™'*). (28) 
/j 

This corresponds to an energy density that has only whatever slow time variation 
the mass-squared-weighted sum of the squares of the absolute values of the tpu's 
may have, and an isotropic pressure that oscillates at the frequencies 2m j and has 
a time average that is zero to this order of approximation (though at the next order 
there are small time-independent pressure gradients that hold up the self-gravitating 
energy density in the approximately periodic or quasi-periodic cases). 
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Directly from the nearly- Newtonian metric ()16|) itself, with the spatial derivatives 
of W being negligibly small, the leading-order linearized Einstein tensor components 
are 

G 00 ~ 2c 2 V 2 U } (29) 
Goi « 0, (30) 
Gij « -2W5 ir (31) 
Thus the Einstein equation in the nearly-Newtonian case becomes 

c 2 V 2 U^Y, m2 i^J^Ji ( 32 ) 
ij 

W*-\ E&he~ 2imit + $je 2imit ). (33) 
8 u 

In summary, for nearly-Newtonian self-gravitating configurations of real self- 
gravitating minimally coupled massive scalar fields, the Einstein-Klein-Gordon equa- 
tions become the coupled approximate partial differential equations (|2Uj) and (|32|h 
which are the time-dependent Newton-Schrodinger or Schrodinger- Newton equations 
[TTJ H2J EH El, plus the additional algebraic equation (|3^|l for the small rapidly 
varying term W in the nearly-Newtonian metric (jlfij) . The conditions for these 
nearly-Newtonian equations to be valid are that J2u |V^j| 2 ^ T 1^1 ^ 1> an d that 
the spatial derivatives of the ipi/s are small in comparison with their typical values 
multiplied by their natural frequencies mj. [Then the Schrodinger equation (}2T))) 
implies that the time derivatives of the ipi/s are also small in comparison with their 
typical values multiplied by their natural frequencies mj.] 

Unless explicitly stated otherwise, in this paper we shall assume that the metric 
is asymptotically flat with asymptotically Lorentzian coordinates (except for a possi- 
ble rescaling). That is, we assume that U goes to a time- and direction- independent 
constant at spatial infinity. Eq. ()32j) then implies that the ipu's must all aymp- 
totically tend to zero at spatial infinity. Eqs. ({20)) and (|3*2*|) are invariant under 
shifting U by a constant or a function purely of t, provided the ipi/s are shifted by 
the appropriate phase factor that is also a function purely of t. (This is simply the 
gauge transformation of replacing the time coordinate t with a new time coordinate 
t' that is purely a function of the old time coordinate t, and of rescaling the spatial 
coordinates appropriately.) One could thus set the asymptotic value of U to be 
zero, making the coordinates asymptotically Lorentzian without any scaling factors. 
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However, in some cases it is more convenient not to make this restriction, such as 
when the time dependence of the tpi/s is purely by a time-dependent phase factor, 
in which case one can cancel the phase factor and make the ipi/s independent of 
time by an appropriate nonzero but time-independent asymptotic value of U. 

For the rest of this section, we shall take the approximate Newton-Schrodinger 
equations (|2()jl and ()32j) as exact and so write = signs rather than rs equal signs. 
However, we must bear in mind that these are actually only approximations, valid 
in the nearly-Newtonian limit, for the actual Einstein-Klein-Gordon equations taken 
as fundamental in this paper. 

The time-dependent Newton-Schrodinger equations (}20j) and (I32j) . now tem- 
porarily re-interpreted as exact equations, may be derived from the classical action 

(cf. ns|) 

^ / dttfxYttrru^u^u - M u ) - c 2 (W) 2 - c 2 |V^/j| 2 - 2m 2 I U\^ IJ \ 2 ] 

dtL NS (34) 

(in units of time squared, to be multiplied by c 5 /hG if one wants a dimensionless 
action), where the classical Lagrangian (with units of time) is 

L N s = I d3x H i™>i$ijl>iJ ~ i/>ij$ij) - E v - E K - E v 

87rc d J 7^ 

j d 3 x imi{$ijipij ~ ^ij^ij) ~ MUoo - E v - E K - 2E P 
J u 

imi{ipijipij — ipijipu) — MUoo — E. (35) 



' NS — 

87TC J 



^7TC 3 



^7TC 3 



IJ 

Here the asymptotic mass of the configuration (essentially the total rest mass, 
uncorrected for gravitational binding energy), in time units, is 



M 47rc3 - ij 
with the rest mass density, in conventional units, being 



y j <PxY. m]Wu? = | / d 3 xp (36) 



p—M^r. or) 

When the Newton-Schrodinger action I^s given by Eq. (}3"4*|) is extremized, so 
that the Newton-Schrodinger equations (|2(J|) and (J32|) are satisfied, one can readily 
see that the asymptotic form of the Newtonian potential is 

U ~ C/oo - — , (38) 
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where M is the asymptotic value of the mass, given by Eq. (|36p above (in time units, 
G/c 3 times the mass M* in conventional mass units, so Mc has units of length and is 
half the Schwarzschild radius corresponding to the mass), and r is the radial distance 
from the center of mass, also in units of length. 

By Eq. (120 j) and the asymptotic boundary conditions given, the mass M is an 
approximately conserved quantity [approximate only to the extent that Eq. ()20j) is 
approximate], corresponding under the approximations used to the exactly conserved 
ADM mass of the spacetime. In some sense it is more nearly the rest mass of the 
matter (in time units), but since only the zeroth-order approximation is being used 
for the rest mass density and for the spatial volume element, M is only a zeroth-order 
approximation for the rest mass as well. 

The various other energies appearing in the Newton-Schrodinger Lagrangian ()35j) 
(smaller than M by factors that in equilibrium are of the order of a typical value of 
\U\, which must be much less than unity for the nearly- Newtonian approximation to 
be valid) are the following quantities (in time units, inserting the factor of 1/c into 
the right hand side of each of the equations with squares of spatial gradients, and 
the factor of 1/c 3 into those without them, under the assumption that the spatial 
distances and gradients are being measured in length units rather than in the time 
units that would avoid the need for all factors of c): the (positive) Newtonian 
potential gradient energy 

Eu = -J— [ d 3 x(VU) 2 , (39) 

57TC J 

the (positive) scalar field gradient energy or matter kinetic energy 

E K = ^- I d 3 xJ2Wu\ 2 , (40) 

the (indefinite in sign, though negative in static equilibrium) matter potential energy 

E v = -± 1 [ d 3 xJ2^ 2 iU\^u\ 2 = -J d 3 x P U, (41) 

the (indefinite in sign, though also negative in static equilibrium) rescaled gravita- 
tional potential energy 

E P = l -E v - l -MU OQ = -^-J 'd 3 xJ2™ 2 i(U -U^ul 2 

Zj oTTC J j j 

= ^Jd'xpiU-U^), (42) 
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and the total Newtonian energy 



E = Eu + Ek + Ev- MUoo = E v + E k + 2E P 

= ^Jd 3 xY / [c\^U) 2 + c 2 \V^ u \ 2 + 2mj(U-U oc )\i: IJ \ 2 ]. (43) 

The extrema of the Newton-Schrodinger action I^s given by Eq. (J34j) . with U 
fixed to a constant value at spatial infinity, and with the ipjj functions falling 
off sufficiently rapidly at spatial infinity, are solutions of the Newton-Schrodinger 
equations ({H and l|32|) . 

When these equations of motion are satisfied, one can readily show (by integra- 
tion by parts, etc., using the boundary conditions given in the previous sentence) 
that M is a constant of motion (as mentioned above), that 

E v = MU^ - 2Eu, (44) 

that 

E P = -Eu, (45) 

that 

E = E K + E P = E K - E v , (46) 

and that the various energies (all of which have absorbed a factor of G/c 3 to have 
units of time) have the following time derivatives: 

E P = X -E Y = -E v = -E K = ^J d 3 x J ■ VU, (47) 

where 

J = Eh^fcV^ - V*jV&j)] (48) 

is the mass-current flux vector in the conventional units of mass per area per time 
(when we remember that m; has units of inverse time and that the ipu's are dimen- 
sionless) . 

One can then see from Eqs. (J46|) and (|47|) that the Newton-Schrodinger equations 
(J2UJ) and (J32J), along with the asymptotic boundary conditions, imply that the total 
nonrelativistic Newtonian energy E is conserved, or constant in time. Using the 
expressions flBUj) and (|4~0Jl above, and using Eq. (jJBJ), we can see that Eq. (|4^j) may 
be written, for E in time units (or E* = c 5 E/G in conventional energy units), as 



c" 



E = ^ = ^ / d 3 x (Y, Wij\ 2 ~ (Vf/) 2 ) = const. (49) 
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This is a first-order correction to the rest mass energy of the configuration in the 
total ADM mass-energy. [But since Eq. (|3H|) for M is only correct for the rest mass 
to zeroth order, the correct first-order expression for the ADM mass is not simply 
M + E.} 

If one has a static solution of the Newton-Schrodinger equations (jSDj) and (jH2j) . in 
which U and each ijjfj is constant in time, then these are extrema of Ejj+Ex+Ey = 
E+MUoa. To put it another way, they are extrema of the total Newtonian energy E 
with the constraint of fixed M. If one finds these extrema by the method of Lagrange 
multipliers, extremizing E — AM, then the Lagrange multiplier is A = — Uoo, which 
thus is dE/dM for a continuous sequence of such extrema. 

If one takes such an extremum static spatial configuration and replaces U(x) by 
U(x) = U(e a x) and each ipu(x.) with ipu(x) = e^jj^x) for small constants a 
and (3, then for the variation of Ejj + Ek + Ey — E + MU^ to vanish to first order 
in both a and (3, and with the use also of Eqs. (JUJ-fjlBJ), one can readily see that 
the static equilibrium configurations have 

2 

E v = 2E K = -2E V = -E P = -MU^ = -2E. (50) 

The relation that the gravitational potential energy is twice the negative of the 
kinetic energy, Ep = —2Ek, and hence that E = —Ek, is just the usual virial 
relation for an attractive inverse-square force. The relation that E = —(l/3)MU OQ 
involves the scaling behavior of the Schrodinger equation and implies that for each 
ipjj to be static (e.g., not to have any time-dependent phase factor), one must have 
Uoo > (essentially for the system to be bound). 

A solution of the approximate time-dependent Newton-Schrodinger equations 
(120)1 and ()3*2*|) is determined by an initial spatial configuration of the ipi/s (and, for 
the quantities in those equations, but not for the gauge-invariant physical quantities, 
by the gauge choice of the asymptotic value of U). If the scalar fields all disperse 
indefinitely, then U will tend to U^ everywhere in space asymptotically with time 
so that the integral of the negative — (VU) 2 term in the integral (|49)) goes to zero, 
and thus the Newtonian energy E must be non-negative. 

But if Eq. ()49j) gives a negative Newtonian energy E, then the scalar fields cannot 
all disperse. To the extent that Eqs. (|2*0*|) and (|52*|) are accurate, at least some of the 
scalar field energy must remain gravitationally bound indefinitely. Unless the scalar 
fields collapse gravitationally into configurations violating the nearly- Newtonian ap- 
proximation being used (and perhaps leading to continued gravitational collapse 



16 



into one or more black holes), the scalar fields will continue to oscillate, giving an 
oscillaton. 

Negative Newtonian energy E is thus a sufficient condition for an oscillaton (or 
for gravitational collapse as a possible alternative) in the nearly- Newtonian limit, 
though it is not a necessary condition, as one can have an initial configuration 
which asymptotically tends to a bound part with negative Newtonian energy and a 
dispersing part with a larger positive Newtonian energy. (See [2] for an example of 
this.) 

In any case, we see that, at least in the nearly-Newtonian case, gravitationally 
bound oscillatons are quite generic (unless all but a set of measure zero collapse 
gravitationally into black holes, which naively seems unlikely), occurring if the single 
inequality E < is true (and also in other cases when just part of the system 
is bound). That is, any perturbation of a nearly-Newtonian oscillaton within a 
sufficiently small neighborhood in the space of perturbations gives another oscillaton 
(unless it collapses). In this sense nearly-Newtonian oscillatons are apparently stable 
under small perturbations, to the degree that Eqs. PH|) and are accurate and 
continue to remain valid (e.g., when one ignores the scalar field radiation considered 
below, and when no continued gravitational collapse occurs that takes one outside 
the validity of these equations). 

Of course, if one requires an oscillaton of fixed total mass (which itself is a 
continuous variable classically, though quantized when one includes the fact that 
the scalar field particles are quantized) to be periodic with some definite period (say 
in proper time at spatial infinity, to make the period gauge invariant), then one gets a 
nonlinear eigenvalue problem with presumably only a discrete set of eigensolutions 
for each mass (modulo gauge transformations, including spatial translations and 
rotations and the gauge transformation of shifting U). Thus I would expect there to 
be discrete periodic oscillatons embedded in an open set of nonperiodic oscillatons 
for a given total mass. 

However, if we include the effect of the oscillating term of the nearly-Newtonian 
metric (jl5j) . the W term given by Eq. (jBHl), we shall see in the next Section that 
it leads to classical emission of each scalar field of mass mi at frequencies roughly 
2m x ± mj, so that generic oscillatons are not classically stable but radiate away 
their energy and scalar fields. The only exceptions appear to be field configurations 
that have a U(l) invariance, so that their total stress-energy tensor is independent 
of time and the W term vanishes. 
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5 Classical Emission of Scalar Fields for Generic 
Situations 



The oscillating W term of the nearly-Newtonian metric (JT5j) gives, in the Klein- 
Gordon equation of the scalar field transitions from its dominant frequency 
near mi to frequencies near 2m^ ± mi that radiate away, carrying off energy and 
causing a generic oscillaton to decay Here we assume that the terms given in Eq. 
(J 17)) describe the non-radiating scalar field oscillations of an oscillaton. 

To describe the radiation, extend Eq. (fTT)) to include scalar field oscillations at 



these emitted frequencies, so the dimensionless rescaled scalar field 4>u = y8irG/c 2 <3>/j 
has the form 

<t>u « V/je" <rojt + $/je <mjt 

+ J2(XiJKe- i{2mK+mi)t +XiJKe i(2mK+mi)t + X ^ 

K 

(51) 

where not only the ipi/s and their complex conjugates, but also the Xhk's and 
x'ijk s an d their complex conjugates, are functions varying much more slowly than 
the frequencies mi and Since fields that can radiate away must have frequencies 
larger than mj, in the sum above we can omit the x'ijk s which have mx < mi. 

Of course, the sum over K includes /, so the radiation field occurs even if there 
is only a single mi or even just a single real scalar field (one choice for the indices 
I J denoting the field). (Then we can omit the x'ijk an d x'uk terms, since their 
frequencies would be just mi, and they would just give small corrections to the ipu 
and ipu terms that are not radiating away by assumption.) 

In principle one should include a whole infinite series of frequencies by adding to 
mi all possible positive and negative integer multiples of all m^-'s, but the additional 
terms will generally be smaller yet, and so to a good approximation it is sufficient 
to consider only the ipu, Xuk, and x'ijk terms and their complex conjugates. 

Now the term with the e " ! l 2mjf+m ''' time dependence in the Klein-Gordon equa- 
tion ((TJ for the field <pu in the metric (|TB|) with the oscillating W term included, 
given by Eq. ()33)h yields 

3 

[V 2 + Am K (m K + mi)]xuK ~ --mim K ^ ^ 2 KL ^U, (52) 

4 L 
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and the term with the e l \ 2mK mi ) t time dependence gives 

3 

[V 2 + Am K (m K - mi))x'u K ~ j m i m K ^kl^u- ( 53 ) 

4 L 

We shall assume, except where explicitly discussed below, that there are no scalar 
waves incoming from infinity, so we shall impose outgoing wave boundary conditions 
on the Xuk's and x'ijk s - Then they have the form 

"XrnTmj/ r P 2i \J m K (m K +m7)|x-x'| 

Xuk(x) « / d 3 x'^ _ £ ^MlM*), (54) 

1m r1T) is r p 2iy/m K (m K - mi )\x-x'\ 

x'ijkM - / ^ — E^(x')^(x'). (55) 

These represent scalar waves that are propagating outward at asymptotic speeds 



2Jm K {m K + mi) 

vijk = — 56 

2rriK + mi 

and 

2^m K (m K -m/) 
= 2m K - mi (57) 
respectively, which generically are within a factor of the order of unity of the speed 
of light (taken to be unity in these equations). 

In contrast, the ipi/s are assumed to be localized almost entirely within some 
region that we shall call the system. Since the i/ji/s are assumed to be slowly 
varying with respect to the natural frequencies mi, we assume that the ipi/s within 
the system do not change much during the light travel time across the system, or 
during the time it takes for the waves represented by the Xuk's and x'ijk s to 
traverse the system. That is why we can use an instantaneous approximation for 
the propagators in the formulas (J54*j) and (J55j) for the Xuk's and x'ijk s - 

If we surround the system by a sphere much larger than the dominant region over 
which the ipifs are significant, then we can calculate the flux of mass out through 
that sphere in the scalar waves represented by the Xijk's and x'ijk s gi yen above. 
When we average this over the oscillation periods and do a bit of algebra that is not 
repeated here, we get a classical mass loss rate of the oscillaton that is 

dM . , dmjmj^ r d 3 xd 3 x' 

^ IJKLM 256VT 2 J |X-X'| 
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x (2m K +m I ) sin(2Y / mx(mK+m/)|x-x'|)^ L (x)^/j(x)^ M (x')^/j(x / ) 



+ (2m K -m I ) sin(2^ m K (m K - mj) |x -x' | )ip 2 KL (x)^/j (x)^ M (x')^j j (x') 



(58) 



In the second of the two terms inside the integrand, whenever — m/ is negative, 
that term is to be omitted, since the corresponding x'ijk * s spatially exponentially 
damped rather than having the oscillatory behavior representing an outgoing wave. 
(Of course, if mx — mi = 0, this second term is just zero, so we need not consider 
it in that case either.) 

Remembering that I am using units and conventions in which M has the units 
of time, in which the m/s have the units of temporal frequency (inverse time), in 
which the ipi/s and their complex conjugates are dimensionless, and in which either 
the spatial coordinates x l (represented above by the 3-vector x) have the units of 
time or else the speed of light is set equal to unity, it is easy to see that both sides 
of Eq. (JHHJl are dimensionless. 

Since the classical mass loss rate formula (}58|) is rather complicated when there 
are several real massive scalar fields of different masses, it may help to give it explic- 
itly when there is only one real scalar field of mass (or actually natural frequency) 



In particular, when ipi and ?/> 2 each have the same magnitude everywhere and are 
everywhere 90° out of phase, or ipi = ±i-0 2 (so the two real fields can be interpreted 
as forming a single complex field with a global U(l) symmetry), then the mass loss 
rate is zero. This is a case in which the stress-energy tensor of the scalar fields does 
not have an oscillating component, and so there is no oscillating W term in the 
metric (fTHJl . 

There is also the analogous case in which one has an arbitrary number (greater 
than one) of fields of each mass mj, when one has zero everywhere for the sum of 
the squares of the ^jj's (not the squares of the absolute values of these quantities) 



m: 




(59) 



When there are two scalar fields of the same mass m, then one gets 




x [$(x) + ^ 2 2 (x)p 2 (x') + ^ 2 2 (x')]^i(x)^i(x') + ^ 2 (x)^ 2 (x')]. (60) 
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for each fixed / (for each different mass mi, the sum being over J that labels 
the different fields with fixed I and hence with fixed mass mi). Then again the 
stress-energy tensor and the metric (fT6*j) has no oscillating term, and so there is no 
generation of outgoing Xijk or x'ijk waves. 

In these cases, if the ipifs are not stationary (or stationary up to a time- 
dependent phase factor) but are slowly changing their form, then although there 
may be no oscillations at the frequencies mi or their sums or differences, the metric 
would still have a slow time dependence, and this would presumably lead to some 
scalar field radiation, though at an amount presumably considerably reduced from 
what it would be the case if the sum of the squares of the ipifs for at least one / 
were different from zero so that there would be the more rapid oscillations in the 
metric. 

The only case in which I would expect absolutely no scalar radiation would be the 
case in which the metric is absolutely stationary. Otherwise it would seem extremely 
unlikely that the outgoing radiation at all possible multiples of the metric oscillation 
frequency, plus or minus the natural frequency of the fields that can potentially be 
emitted, would be zero. However, I have not tried to find a rigorous proof that there 
are not exceptional counter-examples to this conjecture. 

6 Classical Emission of Scalar Fields with 
Spherical Symmetry 

A simple subset of the set of all oscillatons is the set in which the metric and all of 
the scalar fields have spherical symmetry. In the nearly-Newtonian limit (which I am 
taking to exclude gravitational waves), the spherical symmetry of the metric follows 
from the spherical symmetry of the scalar fields, and the spherical symmetry of the 
scalar fields (including the outgoing waves) follows from the spherical symmetry of 
the ipifs. Therefore, we basically just need to assume that each ipu = ipu(t,r). 

The Klein-Gordon equation (0) for the scalar field implies that ipu obeys the 
Schrodinger equation (j2*Uj) . which for the spherical symmetric case becomes 

i(? 

ipu ~ 7; (ripu)" - imiUipu, (61) 

Zmir 

where, except for the prime on x'ijk an d cm dummy variables inside integrals, a 
prime henceforth denotes a partial derivative with respect to r (or later, with respect 
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to a rescaled radial variable x = kr/c). One can avoid the explicit r's in this equation 
by defining 

fjj^^^riPu^r), (62) 



which makes the Schrodinger equation take the form 



IC 



2mj 



fij ~ imjUfu. 



(63) 



The Newtonian part of the Einstein equations, Eq. (J32j) . can be integrated in 
the spherically symmetric case to give at each time 



1 pr poo 

err Jo Jr> jj 



(64) 



Here of course the primes and double primes on the r's inside the integrals just 
denote dummy variables to be integrated over, and not derivatives with respect to 
r. 

The spherically symmetric analogues of and for the Xuk's and x'ijk s 
obeying Eqs. (J52J) and (JoHj) (to solve the Klein-Gordon equation) are (with c = 1 
for simplicity here and in many formulas following) 

"SmirriK 

XiJK{r) r > 



8WmK(mK + mi) r 



e 2i^ mK (m K+mi )r j dr y gin (^ mK ( mK + uij) r ) ^KL{r')ipu(r') 



f dr'r' sin (Jm K (m K + mi) (r' - r))J2 ^KL{ r ')^ij( r ') 



(65) 



x'ij K { r ) 



x 



3mimK 



\\JmKijnK — mi) r 

eVmifK-miJr f°° dr > r > sin (^ % ( % _ m/ ) r ) £ ^ {r')$u{r') 

J 7 

(66) 



dr'r' sin {\jm K {m K - mi) (r - r)) ^ ^K L ( r ')' l l'ij( r ') 

L 



The classical mass loss rate becomes 



E 



§m\m\ 
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X 



+ 



2rriK + rrij 



2\JmK{niK + mi 

2rriK — ttij 
2yJm K (m K - m 7 ) 



/ drr sin (2 J m K (m K + mj)r) V ip 2 KL {r)ipu{ 
Jo i 

drr sin (2^m K (m K - mj)r) y £ J ^K L { r )^ij{'' 



(67) 



As in Eq. (|58|) . in the second term that has the factors of ym^m^ — mi), if 
mK — vni is negative, the corresponding term in the sum over modes is to be omitted, 
since it corresponds to a x'uK~ m0< ^ e frequency 2mx — mj that is below the natural 
frequency mj and hence to a mode that is not freely propagating at large radial 
distance to carry off mass. 

When only one real scalar field of mass m is present, the classical mass loss rate 
becomes 



dM 27m 4 f°° , . r . ,„ . 
— « - .„ / drr sin V8mr w (t") 

With two scalar fields of identical masses present, one gets 

dM 27m 4 



(6* 



dt 2 n /2 



drr sin (V8mr)(ipl +^2)^1 



2 




+ 


/ 

.A) 



drr sin (v/Smr)^^^)^ 



(69) 

Again one can readily see that if the sum of the squares of the complex if) fields with 
identical masses are zero, then there is no classical mass loss, at least at this level 
of approximation. With just two fields of identical masses, this condition is that 



7 The Simplest Spherical Oscillatons 

Now let us focus on the simplest case, in which there is a single massive scalar field 
in a finite-mass spherically symmetric nearly-Newtonian configuration that is very 
nearly periodic in time and has no nodes. In particular, require that ip be spherically 
symmetric, independent of time (except for a possible slowly varying phase factor), 
and nowhere zero (no nodes, though asymptotically zero at spatial infinity). 

By a suitable choice of the hypersurfaces of constant time, one can cancel the 
phase factor to make ip real and positive everywhere, which is what I shall assume, 
at the cost of having U approach a nonzero constant at spatial infinity. Then the 
time-dependent Newton-Schrodinger equations (|2Uj) and (J32J) for a single scalar field 
with a single real ip become the time-independent Newton-Schrodinger equations 



23 



[TH IT2*1 IT3*1 IT4*] (where for simplicity I am using units in which c = 1 so that I can 
drop many occurrences of factors of c that one can easily put back in by dimensional 
analysis if needed with other choices of units), 

V 2 V ~ 2m 2 Uip, (70) 

and 

V 2 U^m 2 if) 2 . (71) 

In the case of spherical symmetry which I am now also assuming, the equa- 
tions take the form (with if) real) of two coupled second-order ordinary differential 
equations, 

if;" + -if;' « 2m 2 Uip, (72) 
r 

and 



U" + -U'k m 2 if; 2 . (73) 
r 

As noted previously fTJ EI| , this system of equations has the scale invariance 

(ip, U, r) (-> (X 2 ip, X 2 U, A"V). (74) 

Inserting explicitly the speed of light c so that the radius r can have units of length 
instead of time as it implicitly does above, I shall set 

r = |, (75) 

u - w < 77 > 

so with k having the units of frequency or inverse time, x, S, and V are the dimen- 
sionless variables used by ^Hl (except that what I now call x, they call r). Note that 
this V, which is just a rescaling of the Newtonian potential U with its sign reversed, 
is not to be confused with the V(t, x k ) in the metric ([14)1 . which is negligible in the 
nearly-Newtonian metric. 

Although so far I have used a prime generally to denote a derivative with respect 
to the radius r, when I am using the dimensionless rescaled radial variable x instead 
as the independent variable, I shall use a prime to denote a derivative with respect to 
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x. Then Eqs.()72|) and ()73|) (where the prime denoted d/dr) become the dimensionless 
time-independent Newton-Schrodinger equations 



( X S)" = -xSV (78) 

and 

(xV) n = -xS 2 . (79) 

Note that I have now replaced the ~ signs with = signs, even though these equations 
are only an approximation to the actual Einstein-Klein-Gordon equations, though 
an approximation that becomes arbitrarily good in the nearly- Newtonian limit. 

These time-independent Newton-Schrodinger equations are the same as Eqs. (6a) 
and (6b) of except for the replacement of their radius r by my dimensionless 
rescaled x. The fact that the arbitrary constant k does not appear in these equations 
illustrates their scale invariance. 

Another form of the equations that is helpful for some of the analysis below is 
to use 

p=l(S-V) (80) 
instead of V in Eqs. (f7%|) and (J7HJ), which then become 

(xS)" = -xS{S - 2p) (81) 

and 

(xp)" = xSp. (82) 

Yet another set of variables to use that were particularly convenient for numerical 
integration of the equations and for the interpretation of the results are 

X = xS, (83) 

M = -x 2 V, (84) 
M 

w = -(xV)' = V. (85) 

x 

The variable X is the same as that used in j!4j . The asymptotic value of M.(x), 
say M.OO-, is what is called B in ^3] arid what is called / in Eq. (2.19) of [Hj. The 
asymptotic value of w(x), say w^, which is positive, is by Eq. (J83j) the same as the 
asymptotic value of —V(x), say — V^, which is the same as —A in [15] . 

Ai(x) has the interpretation of the rescaled mass interior to the sphere at x. 
From the fact that asymptotically one has U ~ — M/r and V ~ + Ai/x 
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« = ^ = ^. <■») 



(assuming a finite-mass oscillaton in which the mass-energy density, proportional to 
S 2 , asymptotically rapidly approaches zero), one can see that the mass (in units of 
time) interior to a sphere of radius r = cx/k is 

M(r) = ^M(x). (86) 
The variable w can be written as 

2m 2 

w = —{U + rg), (87) 

where g is the acceleration of gravity at the same radius r where the gravitational 
potential U is evaluated, 

dU _ M{r) 
dr r 2 

Thus w can be interpreted as a rescaled value that the gravitational potential U 
would have at twice the radius r if U had a uniform gradient from r to 2r, with 
this uniform gradient having the same value of the actual gradient at radius r. If 
the mass-energy density dropped precisely to zero outside some radius, w would 
be constant outside this radius, at the value w = —Voo = (2m 2 /A; 2 )f/ 0O . In the 
actual case in which the mass-energy density, proportional to S 2 , drops exponentially 
toward zero, w is exponentially close to — and is therefore a good integration 
variable to use to evaluate the asymptotic value of V, that is V^. 
The differential equations in terms of these variables are 

M 

X" = ( w )X = -VX, (89) 

x 

M' = X 2 , (90) 

/ X 2 , , 

w' = — . 91 

x 

The initial conditions for these variables (initial in x, of course, not in time, 
since all of the quantities being considered presently are independent of time in the 
approximation that we initially ignore the mass decay rate) are that at x = 0, we 
have X = 0, X' = S , M = 0, and w = -V . 

As discussed in Ej , for an everywhere regular static spherically symmetric 
solution to the Newton-Schrodinger equations, S and V must be smooth everywhere, 
and S must be decreasing exponentially at spatial infinity (x — > oo). At the origin 
{x = 0), S and V must have finite values, Sq and Vq respectively, and must have 
zero dimensionless radial derivatives, S' = and V = 0. 



26 



Because of the scale invariance, the only independent nontrivial parameter for a 
solution regular at the origin is the ratio Sq/Vq. Integrating out from the origin gives 
a solution that diverges at finite radius (with S going to +00 and with V going to 
— 00 there) if Vo < 0, so we shall choose V Q > 0. Using the scale invariance, without 
loss of generality we can and shall set Vo = 1, leaving the nontrivial parameter to 
be Sq. If So = 0, then we just get the trivial solution S = 0, V = Vq, which is 
flat spacetime with no matter, a solution we shall discard as previously studied by 
other people. By the symmetry of the equations under S 1— > —S, we can thus choose 
S > 0. 

If So is too large, Eq. (|T9|) implies that xV (initially growing as x) bends down 
rapidly, so that V goes negative while X = xS is still growing. (Initially X = xS 
also grows linearly with x, as Sqx, but, like xV, it also bends down. However, Eq. 
(J7HJ) implies that it does not bend down so fast as xV bends down, for S > V > 0.) 
Then when V becomes negative, X' grows with x, and so X grows faster and faster, 
and eventually so does S, by Eq. (|75j) or the equivalent Eq. (jSHj). Eq. (|7§|) imples 
that then xV and eventually also V gets more and more negative. In fact, then 
S goes to +00 and V goes to —00 at a singularity of infinite mass at finite x. Of 
course, we shall discard these solutions. 

On the other hand, if Sq is positive but too small (e.g., less than Vq = 1 [TI]). 
Eq. ()79j) implies that V will stay positive long enough for Eq. (|78jl to imply that 
x S will oscillate (with characteristic period in the x- variable of 2ir/^/V if V were 
constant). However, we want the regular solution with no nodes, the solution with 
the largest value of S that does not lead to a singularity. This value of 5*0 is an 
eigenvalue for the system. 

At this eigenvalue, X = xS will bend over from increasing at x = (X' = Sq 
there) to decreasing again toward zero value asymptotically, but never crossing zero. 
xV will also bend over from increasing (at unit rate) at x = (since we have chosen 
Vo — 1) and will cross zero to become negative and will eventually keep decreasing at 
the asymptotically rate given by the asymptotic negative value of V, say V^. That 
is, S will start out at So with zero derivative with respect to x (S' = initially) but 
will then bend down to approach zero asymptotically at large x (while then bending 
back upward just enough to keep from crossing zero, but never quite leveling out, 
except asymptotically). Similarly, as a function of x from x = 0toa: = oo,V will 
start out at Vo = 1 also with zero derivative with respect to x (V = initially) but 
will then bend down and cross zero before bending back up to level off toward some 
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negative constant asymptotic value V^. 

For this eigensolution, since X starts at zero and initially increases linearly with 
x, and since M.' = X 2 with Ai(0) = 0, and since w starts at — Vq = — 1 with 
w' = X 2 /x, we have initially (near the origin, a; < 1), 

M ~ \sy - ±S 2 x 5 , (92) 

w ~ -1 + ^S V, (93) 

X ~ S x - ls x 3 , (94) 
o 

S~S - ^S x 2 , (95) 

F ~ 1 - l -S 2 x 2 , (96) 

p~^ -l)(l + ^ x 2 ). (97) 

Then as x is increased to oo, X will at first increase, while bending downward 
and eventually passing a maximum and then decreasing. While decreasing, X will 
pass an inflection point (at the point at which V crosses below zero) and will then 
bend upward to level out asymptotically as it also approaches zero asymptotically. 
At the same time (here meaning during the same evolution in x from to oo), M. will 
start from with zero slope and curvature and initially grow as the cube of x (i.e., 
as the volume interior to the sphere of radius r = cx/k) but eventually will reach an 
inflection point (at the point at which X reaches its maximum) and then gradually 
level off to approach its asymptotic value Ai^. The variable w will start at —1 at 
x = with zero slope and will bend up to cross zero, before bending downward to 
level off asymptotically and approach its asymptotic value = —V^ > 0. 

For use below in calculating the quantum decay rate of this type of oscillaton, it 
is also of interest to integrate the variable D(x) given by the differential equation 



X 
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D' = x 2 S* = -^ (98) 
x 2 

with the boundary condition D(0) =0. Since the mass density in conventional units 
is (with m and k in frequency units) 

P 4ttG SyrGm 2 ' 1 ' 
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the asymptotic value of D, namely D^, is proportional to the integrated square of 
the mass-energy density (and hence to the total annihilation rate of two scalarons 
into two gravitons): 



I have used the differential equation routine of Maple 8 to evaluate to high 
accuracy the eigenvalue So and the asymptotic values M.^ and Woo, making many 
overnight and several-day runs on a Dell 2350 PC with the Maple "Digits" parameter 
set to values as high as 62 (but more often 40 or 50) and the "relerr" and "abserr" 
parameters set as small as 10~ 35 . The calculations performed roughly 10 8 function 
evaluations in a one-day run, so several billion function evaluations were performed. 

I did face the annoyance that after completing a calculation, Maple 8 would 
often crash when I attempted to save the file and would erase the entire file (even 
if I had saved the pre- run version). I had to get used to renaming and saving files 
before making runs and to copying down key results by hand, and then copying and 
pasting small pieces of the files before attempting to save the entire file, so that 
I would have something left in case the entire file crashed and disappeared. This 
meant that by the end of the calculations, I had thousands of digits copied down by 
hand, and hundreds of different files saved. 

The differential equation routine handled the regular singular point of the dif- 
ferential equations above at x = 0, so I could simply set the initial conditions there. 
However, it could not handle the essential singularity at x = oo, so I could only 
do the integration out to a finite value of x. Any errors in Sq and in the finite- 
differencing technique would lead to an S(x) that did not continue to approach zero 
indefinitely but rather would diverge away from zero and eventually go to +oo (say 
if So were set too large) or to — oo (say if S were set too small). Figure 1 of [13J 
illustrates this phenomenon. 

Therefore, for large values of the radial variable x, I developed an asymptotic 
series solution in terms of decaying and growing exponentials and in terms of an 
expansion in terms of 1/xf, where Xf is the final value integrated to (generally 
in the range 30-50), for estimating the exponentially-growing error terms (from an 
inaccurate Sq or from numerical errors at smaller x) , for getting rid of their dominant 
effects, and for estimating what the corrected solution would give (through five 
orders in 1/xf) for M.^, w^, and D^. 

Once I had found So to a fairly high degree of accuracy, I assumed that the 




(100) 
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error in X that I could thus estimate at some Xf was linear in the error of So, an 
assumption that fit the data very well when I tried four slightly different values of So 
(e.g., differing at the 35th digit). Then I could interpolate to a value of So that would 
be predicted to give no error. The interpolations to four predicted values of So, from 
four of the six pairs of four trial integrations at the highest accuracy settings, agreed 
to 45 digits, and I shall give my best estimate for So to 42 digits, using a well-known 
answer from elsewhere in choosing this number of digits. (The high agreement of 
the interpolated values for So suggested that the random numerical errors from the 
numerical differencing were much smaller than the values of 10~ 35 I set in these 
multi-day computer calculations for the numerical routine parameters "relerr" and 
"abserr.") 

However, this agreement does not take into account a possible systematic error 
from the numerical differencing. From other trials at larger values of "relerr" and 
"abserr," I crudely estimated that the resulting interpolation for So might have an 
error around 10~ 35 , though it might be smaller. Therefore, although I am giving 
my estimate of So to 42 digits below and believe that value of So is the best value 
(to that number of digits) to use in the routine I used in Maple 8 with the accuracy 
parameters as I set them, I am not very confident about the accuracy of the last few 
digits given if one were able to integrate the differential equations without numerical 
errors. Therefore, the last few digits should be taken as merely conjectural. 

Then from the estimates from each of the four most accurate trials (at very 
slightly different trial values of S near the true value) for M^, w^, and (each 
of which attempted to correct for the exponentially growing error term, which was 
a very small fraction of the total for X at the values Xf — 35, Xf — 40, Xf — 45, 
and Xf = 50 that I used, at least for the trial values of So that were sufficiently 
close to my final estimate), I could also do linear interpolations to what their values 
would be expected to be at the correct value of So (again neglecting a possible 
systematic error from the numerical differencing). The four interpolations agreed 
for .Moo through 36 digits, for through 38 digits, and for through 44 digits 
(more digits in this case because the integrand for D falls exponentially twice as 
fast as the integrands for M.00 and u>oo)- However, because of a possible systematic 
error from the numerical differencing, these interpolated asymptotic values are not 
likely to be accurate to these numbers of digits. Based on a comparison with trials 
at larger values of "relerr" and "abserr," I suspect that they are accurate through 
about 32 digits. Nevertheless, I shall give them to 36 digits, again with the last few 
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digits being just my best conjectures from the long calculations. 

The number of runs at small values of "relerr" and "abserr" was limited by the 
large amount of running time for the Maple 8 program to achieve its acuracy criteria. 
For example, one of the final four runs I did, all four with "Digits" set to 50 and 
"relerr" and "abserr" each set to 1CT 35 , took 3.8 days to get just to Xf = 35 (where 
I terminated the calculation), though I had asked it to give the results at a sequence 
of x/s out to 50. (A final run, using a value of So that agreed to 40 digits with 
my best estimate after using the results of that run, did get to Xf = 50 within 1.6 
days. I'm not sure why it went further and faster.) So on my PC it would be too 
time-consuming for me to try to give the results to much more accuracy than what 
I am giving below. 

The values I obtained were 

S « 1.088 637 079 429 044 995 770 660 314 259 096 461 243 20, (101) 

Moo ~ 3.618 701 237 823 656 810 212 940 073 900 601 34, (102) 
^ « 1-065 731 278 365 451 058 714 678 859 310 396 16, (103) 
« 1.320 680 334 028 957 063 595 721 629 052 089 05. (104) 

As noted above, I think these numbers are likely to be correct to at least 30- 
35 digits, but I am not really sure about the last several digits, which should be 
interpreted as just my best current guess for what they are. 

One can see that my value for So confirms all but the last of the 15 digits given 
for this quantity by [T^]. (And if So were larger by just 4.23 x 10 -18 than it seems 
to be, the last digit of the best 15-digit approximation would have indeed rounded 
up to the 5 that ^3] gave. However, modulo some highly unexpected conspiracy in 
numerical errors, I am fairly confident that my value is correct to at least 30 digits, 
since I could get agreement to more digits even when I changed from numerically 
integrating one form of the set of differential equations to numerically integrating 
another form.) 

From these numbers, one can of course construct various combinations of them, 
such as the scale-invariant quantity A/B 2 that (El discuss: 

A = -- ~ -0.081384 603 921072 995 209 339465 504462 035 7. (105) 
B 2 M 2 ^ 

The first two nonzero digits of this quantity seems to agree with the value ^H] 
plotted in their Figure 4, but they do not list its numerical value. 
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From the scaling relation given by Eq. (j8o|) . we can get the small dimensionless 
mass parameter of the nearly-Newtonian oscillaton, 

li = Mm = ^Moo- (106) 
2m 

Since in the end we want to express other properties of the oscillaton in terms of 
/i, we shall actually invert this to get the scaling parameter k (which has units of 
frequency, as does m) as 

For example, the value of the Newtonian potential at infinity is 

k 2 2w 

= -Voo = -r^V « 0.162 769 207 842 145 990 418 678 931 008 924 071 u 2 

2m 2 Mio 

(108) 

We can then express the fractional binding energy of a nearly-Newtonian oscil- 
laton as 



1 OC 



with 



e = « 0.054 256402 614 048 663 472 892 977002 974690 5. (110) 

3yvl, 



:3C 



If we take the number of scalar particles to be 

N _ M* _ c 5 M 
TiGm ' 

so the conventional rest mass of the oscillaton is iV times the conventional rest mass 
of the scalar field quantum, M* = Nm*, then the conventional total energy of the 
oscillaton is, to first order in /j, 2 , 

E tot = Nm^il + E/M) = Nm*c 2 - N 3 mleG 2 /h 2 . (112) 

This would be the same value for a boson star with a complex scalar field having a 
U(l) symmetry. For that problem the value was calculated by Ruffini and Bonazzola 
[To] over 34 years ago, getting e = 0.1626. Except for the last digit, this corresponds 
to three times the value above. 

One can also compare my numerical results with the 5-place results for a boson 
star by Friedberg, Lee, and Pang jillj (for n = nodes). In terms of my calculated 
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parameters and my numerical results, their calculated parameters would be 

7o = -Vo/Sq = 1/So » 0.918 579 771 804 638 252 351 006 743 940 091 299 098 45, 

(113) 

7oo = w^/So « 0.978 959 194 486 001 441 125 618 083 933 781 532, (114) 

7i = Moo/Sl' 2 « 3.468 256 171 397 572 159 904 340 931 700 591 17. (115) 

They got % = —0.91858, 7^ = 0.97896, and 71 = 3.46826, in perfect agreement 
with my results rounded to five digits after the decimal point. 



8 Classical Emission from the Simplest Spherical 
Oscillatons 

Now that the nodeless spherically symmetric nearly-Newtonian configurations have 
been found (determined by the single scaling parameter \i = Mm), we need to use 
Eq. (|5Sjl to evaluate the classical mass loss rate. By using Eqs. (|75)) . (|7S)l. and 
(j!07|) . we can write the classical mass loss rate, or power emitted in classical scalar 
radiation (dimensionless when the mass M is in time units, i.e., M = GM*/c 2 in 
terms of the mass M* in conventional mass units), as 

27u 8 

P c = -M classical « f F 2 , (116) 

dxx S 3 (x) sin (ax), (H7) 

a = . (118) 

Because the nearly-Newtonian configurations have /j< 1, the /^-dependent con- 
stant a is very large, a ^> 1. Therefore, the sin (ax) factor in the integral ()117|) for 
F oscillates very rapidly and nearly washes out the integral for large a, causing F 
to be very small. 

We can estimate the value of F by the following method of contour integration: 
Since S(x) is an even function of x, the integrand is an even function of x, and so 
the integral along the real axis from to 00 may be replaced by half of the integral 
along the real axis from —00 to +00. Then the real variable x may be extended to 
the complex variable function of which the integrand is analytic except at 

poles of S(z). One may split sin (az) = 0.5ie~ iaz — 0.5ie mz into the first exponential, 
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which drops exponentially along the negative imaginary axis for the complex z, and 
the second exponential, which drops exponentially along the positive imaginary axis 
for the complex z. By splitting up the integral into the corresponding two pieces, 
the first piece may be replaced by a contour integral making a clockwise loop around 
the lower half plane (Im z < 0), and the second piece may be replaced by a contour 
integral making a counterclockwise loop around the upper half plane (Im z > 0), 
which gives an equal contribution. 

Now one of these contour integrals, say the second one (since they each give 
equal contributions), may be evaluated by finding the residues at each of the poles 
of the integrand. There are a series of poles of S(z) running up the imaginary z 
axis. Because of the e mz factor from sin (az) , the dominant residue will come from 
the pole closest to the real azis, say at z = iyo- 

From an analysis fHJ E| of the time- independent Newton-Schrodinger equations 
in their dimensionless form (J78J) and (|79|). one can see that the solutions are analytic 
over the complex z-plane, except for movable (moving if So were changed) double 
poles with coefficients —6. In particular, near the pole at z — iyo, S has the 
asymptotic form 

s(2) ~(db- (119) 

Because it is S 3 (z) that appears in the integral (|117|) for F, which thus has a 
6th-order pole at z — iyo, one must integrate this factor by parts five times, giving 
the 5th derivative of the e mz factor, and hence 5 powers of a, in the dominant term of 
the result. (I shall drop terms with lower powers of a that arise from differentiating 
the z factor rather than the e mz factor, since they involve higher powers of the small 
quantity //.) As a result, one finds that 

F « -1.87ia 5 y e- ayo . (120) 

Thus there is one more parameter that must be determined numerically before 
we can evaluate explicitly the dominant term (for /j, <C 1) in the classical mass loss 
rate, namely the value yo that locates the pole in S(z) at z — iyo that is above the 
real axis but nearest to it. 

To find yo, one can integrate the time-independent Newton-Schrodinger equa- 
tions (J78|) and (J79|) up the imaginary z-axis (after replacing the real radial variable 
x with the complex radial variable z). That is, one can set z = iy and integrate 
these equations from y = to the singularity at y = yo. If one now wants a prime 
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to denote a derivative with respect to y rather than with respect to x as it does in 
Eqs. (|78|) and (|79jl . one merely needs to reverse the sign on the right hand sides of 
these equations. To have only real quantities in the resulting equations, one could 
also replace each x that occurs once on each side of each of these equations with y, 
so that one has purely real equations to integrate along the y-axis. Because both 
S(z) and V(z) are even functions of their argument, they both remain real when 
one evaluates them along the imaginary z-axis. 

An alternative way to do this, getting a single real equation for both real z and 
purely imaginary z, would be to define v = —z 2 to use instead of x or y as the 
independent variable. Then to get the nearly-Newtonian oscillaton configuration for 
real radius, one would integrate from v — along the negative real f-axis, whereas 
to find the first pole at the imaginary value of the radius, one would integrate from 
v — along the positive real v-axis. The dependent variables S(v) and V(v) would 
both be real along both the negative and positive branches of the real f-axis, though 
they are not even functions of v as they are of z. 

Of course, one cannot integrate numerically precisely to the pole, so one must 
instead integrate to a point near it and then use the behavior there to estimate 
where the pole is. To minimize the numerical error in this process, I thought it 
would be better to redefine the dependent variables so that most of them (all but r 
below, where this dependent variable r is not to be confused with the radial distance 
r that is imaginary for this part of the calculation) do not blow up at the movable 
singularity. I used v = —z 2 as the independent variable and the following four 
dependent variables: 

P=\(S-V), (121) 

q = \Z|, (122) 
S' - V 

r = — , (123) 

Then the time-independent Newton-Schrodinger equations ()78)1 and (|79|) (with 
x replaced by z and a prime then denoting d/dz) give the following four coupled 
first-order differential equations: 

I— ■ <-> 
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dq 
dv 



(126) 




(127) 



(128) 



The initial conditions for these equations are that at v — 0, one has 



p(0) 



2 



(129) 



?(0) = /J' 



(130) 



r(0) 



So(Sq — 1) 



(131) 



12 



s(0) = 



(132) 



These dependent variables were chosen so that g goes to zero at the pole, linearly 
in v, and so that p and r have a negligible effect on the evolutions of q and its negative 
derivative s near the pole. 

By defining the variable u = q/y/v, which like q goes to zero at the pole, and 
developing series solutions in terms of u at points sufficiently near the pole that u 
is very small, one can find that the location of the pole at v = y% is given to high 
accuracy (when u <C 1) by 

y « y/v[\ + u + 0.2m 2 - 0.02m 3 + 0m 4 - 0.0018m 5 - 0.00256m 6 + 0{u 7 )}. (133) 

Therefore, I continued to use the Maple 8 differential equation routine to inte- 
grate Eqs. (|125|) - ([128p from v — 0, where I used the initial conditions (|129|) - (jl32j) . 
out to various values very near v — y$ where u was quite small, and then I used the 
values of u and v there to estimate yo, the location of the pole, by Eq. ()133j) . 

Actually, unlike the case of Eqs. ()89j l -(|9i p . which Maple 8 handled with the initial 
values given at the regular singular point x — 0, Maple 8 gave me trouble trying 
to start integrating Eqs. ()125|1 - (jl28j) from their regular singular point at v = 0, so 
I had to develop series solution near that point and start the numerical integration 
from a slightly nonzero value of v, say V\. I needed to choose V\ small enough that 
the series solution was sufficiently accurate at this starting point, and yet not too 
small, or else the effects of the regular singular point increased the numerical error 
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in integrating from that tiny value of v to a more reasonable value. I found that the 
precise value of v\ had very little effect on the first 36 digits of the result for yo if it 
were around v\ = 10~ 16 , so I generally used that starting value for the independent 
variable v. 

I also had to choose a final value of v, say v%, very near the singularity, so that 
Eq. (j!33j) would be accurate in estimating the precise location of the singularity, and 
yet not so near that the numerical errors from integrating too near the singularity 
would be significant. I tried about 20 values for v 2 and chose the value for so that 
the resulting yo did not not vary much with shifting V2 a small amount either way 
from the value I thus chose for it. Again the result for yo showed stability to about 
36 digits when t> 2 was chosen in a suitable range, such as between 14.843 68426 and 
14.843 68427. 

Of course, there is likely to be some systematic error beyond the random error 
that caused the results to fluctuate in roughly the 37th digit when I changed V\ 
and t>2 slightly, so I am not confident that I really have determined y§ to 36-digit 
accuracy, but since I got a preferred estimate to this many digits from my numerical 
calculations, let me give the pole location y Q , like .Moo, w^, and Doo, to 36 digits: 

y « 3.852 750 221 596 529 691 939 909 038 965 636 80. (134) 

Let me repeat my warning that the last few digits of this result are uncertain. 

Now that we have the last parameter, y , that we need to evaluate F by Eq. 
(I120|) . we can go back and plug the result into Eqs. ()116J) and ()118J) to give the 
classical mass loss rate as 

Pc = -Mclassical ~ ^ (135) 

where 

a = v^-Mo^o ~ 39.433 795 197 160 163 093 719 521 465 426 369 6 (136) 
and where 

o3/2o7 

C = — =— ^a 2 w 3 797 437.776 333 014 909 099 808 643 511 322 13 (137) 

As with the other numerical quantities given above, the last few of the 36 digits 
given for a and C are actually uncertain. 
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I should emphasize that this is what I believe to be just the dominant term in the 
classical emission of scalar waves from a nearly-Newtonian oscillaton when /i C 1. 
I would expect this expression to have a relative error of the order of \x. 

To put it another way, one can presumably write 



so that the calculations performed here have given the three leading terms in an 
expansion for In (1/P C ). 

It is beyond the scope of this paper to do the nonlinear gravitational calculations 
to find the classical mass loss rate when [i is not small, but if one neglects the 
corrections, one can make a very crude estimate for the mass loss rate even up to 
the maximum value of /i, say /i max . 

For what are generally called boson stars (stationary spherical configurations of 
a complex massive scalar field whose phase rotates in a circle in the complex plane, 
equivalent to two real scalar fields oscillating 90° out of phase), the maximum value 
of fj, is 0.633 [TT]. For oscillatons of a real massive scalar field, the initial calculations 
P gave /i max ~ 0.6. Urena-Lopez [lj has given /i max — 0.5522, but Alcubierre et al 
[Hj have given /i max = 0.607, so I do not know what the correct answer is, and I have 
not attempted to calculate it myself. However, this uncertainty is no doubt small 
in comparison with the error of applying Eq. (jl35J) or (|138|) to ji = /i max . 

For example, if the coefficient of the O(fi) term of Eq. ()138j) were the same 
magnitude (but of uncertain sign) as the coefficient of the 1//J, term, namely a ~ 
39.4338, then just taking this single term with, say fi = 0.633, would change In (1/-P C ) 
by roughly ±25.0, or a total range of roughly 50.0 for this quantity, giving an 
uncertainty in the mass decay rate by a factor of about e 50 ~ 5 x 10 21 . One might 
hope that the uncertainty is a lot less, but without calculating the 0(/i) and higher 
terms, I don't see how one can be sure. 

Despite this proviso, if we did naively insert the boson star maximum mass 
parameter /z max ~ 0.633 into Eq. ()135|) . we would get a mass loss rate (dimensionless, 
since our masses denoted by M have the factor of G/c 3 inserted to give them the 
dimension of time) of about 8 x 10~ 21 . However, if it could be larger or smaller by a 
factor of roughly e 25,0 ~ 7 x 10 10 , the dimensionless mass loss rate could be as large 




- -21n- -lnC + 0(/i) 



H n 

39.433 795197160163 093 719 521465426 369 6 , 1 

2 In - 

H fi 

-15.149 837 127 888 728 198 841 232 066 408 315 377 + 0(/x), (138) 
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as roughly 6 x 10 -10 or as small as roughly 10~ 31 . 

Similarly, if /i max ~ 0.607, one would get a crude estimates of the mass loss rate 
as 6 x 10~ 22 , and a crude guess that it might be as large as 1.6 x 10~ n or as small 
as 2.5 x 10~ 31 ; if /i max ~ 0.5522, one would get a crude estimates of the mass loss 
rate as 1.2 x 10~ 24 , and a crude guess that it might be as large as 3 x 10~ 15 or as 
small as 4 x 10 -34 . 

In any case, unless the coefficient of /x. in the 0(/i) term of Eq. ()138|) (or actually 
this entire correction term divided by /x) were negative and had a larger magnitude 
than the coefficient of the l//i term, it seems that the dimensionless mass loss rate 
is always less than about 10~ 9 /x. This means that during a one-radian change in the 
phase of the scalar field oscillation (a time t = 1/m), the oscillaton would have lost 
less than one-billionth of its mass. If instead the correction to applying Eq. ()135)) 
to /i = /x max were negligible, then during the 27r-longer period of a full scalar field 
oscillation, even a maximum-mass oscillaton would have lost less than a billionth of 
a billionth of its mass. 

This result shows why the numerical analyses to date have not shown any insta- 
bility of the oscillaton, since its mass loss rate is so low. 

Figure 6 of [S] for /i = 0.5726 shows an apparent numerical mass loss rate of 
about 3 x 10 -9 ~ 5 x 10 _9 /x in dimensionless units, but since this is several times 
higher than the crude guess above for the upper limit for the mass loss rate, and is 
~ 2 x 10 14 ~ e 33 ~ e 58M times larger than what my rashly-interpolated formula would 
predict for that /i, namely P c ~ 1.4 x 10~ 23 , I suspect that the authors are indeed 
correct in attributing this to "a small amount of numerical dissipation still present 
in our numerical method." However, I can't completely rule out the possibility that 
the 0(/x) term in Eq. ()138|) is roughly — 58/i at /i = 0.5726, so that there conceivably 
might be mass loss comparable to that given by Figure 6 of jHJ. 

Later we shall find that for sufficiently small /x, the classical mass loss rate (j!35)l is 
dominated by a quantum mass loss rate. However, when the classical mass loss rate 
dominates, and when the time t 2 —t 1 to evolve from /x x to \ii is approximately 

*2 - h « -^-e a ^ 2 - -^-e a ^K (139) 
aCm aCm 

For example, we could define i(/x) to be the time to decay from /x = /x max to some 
smaller value of /x. Then if /x max — /x ^> ~ 0.00934, then the magnitude of 

the first term on the right hand side of Eq. ()139|) (with /X2 = /x) is much greater 
than the magnitude of the second term (say with /xi = /x max ~ 0.607). Since this 
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is necessarily the case for /i < 1 where Eq. (|139|) is applicable, we may then drop 
the second term and say that the time to decay down to fi <C 1 from /i max by the 
classical emission of scalar radiation is 

« -^e a/ ». (140) 

9 Precisely- Periodic but Infinite- Mass Oscillatons 

Although it is beyond the scope of the present paper, which just gives numerical 
results for /x <C 1, it would be of interest to be able to calculate the function 
for all (jl < /i max . To calculate 'exact' results (i.e., exact up to numerical errors in 
solving the differential equations), one needs a precise definition of /i max and of /x(t) 
for the time thereafter. 

It is rather hard to define /i max precisely (and the initial oscillaton configuration 
that gives this maximum mass), since any initial configuration is losing mass (as- 
suming boundary conditions of no incoming scalar waves), so one could start with 
a wide variety of initial configurations. But essentially one would like to start with 
one out of a set of initial configurations that lose mass as slowly as possible for each 
initial mass, and then choose the maximum-mass element from the set members that 
do not have rapid mass loss (e.g., at a rate roughly given by the dynamical timescale 
1/m). If the initial decay timescale is of the order of [/i^ ax /(aCm)]e a/ ^ max , as the 
weak-field formulas would suggest, then to the degree that this is much larger than 
1/m, one can define the maximum-mass initial configuration to that accuracy, i.e., 
with a relative error that would be expected to be of the order of (aC//i^ ax )e~ a//Mmax . 
However, to define it to more accuracy (e.g., 'precisely') would require more care. 

Although I do not see any highly preferred way to make the definition precise, I 
can propose the following ad hoc method, which leads to a consideration of precisely- 
periodic but infinite-mass oscillatons: 

Temporarily relax the condition of no incoming waves that has been fundamen- 
tal to the discussion so far. Then it appears that one can have precisely periodic 
oscillatons, though now with infinite mass, so that the spacetime is not quite asymp- 
totically flat (although having curvatures falling off fairly rapidly, so, for example, 
the spatial integral of the Kretschman invariant, the total 4- fold contraction of the 
square of the Riemann tensor, is finite at any time). 

The idea of these precisely periodic (though infinite mass) oscillatons is that 
not only does the oscillation of the scalar field at frequency m (with respect to a 



40 



suitably scaled coordinate time t in a gauge or choice of coordinates in which = 
and goo is independent of t) have precisely the right phase to drop exponentially to 
zero at large spatial distance by the gravitational binding of that mode, but also 
the oscillations at all odd multiples of m also have precisely the right phases to 
drop exponentially to zero at large spatial distance by the gravitational binding of 
those modes as well, when one takes into account the gravitational field not only of 
the mode at frequency m but also the gravitational field of all the higher-frequency 
modes. 

That is, one has a sequence of relations for the modes at each frequency that 
start as follows for the lowest mode: The mode with frequency m is chosen to have 
the right phase (e.g., the right relation of its initial-in-r value at r = 0, determined 
by the quantity So above, to the the value of goo there, determined by the quantity 
Vo above in the weak- field or nearly-Newtonian limit), so that once the stress-energy 
tensor of the scalar field (predominantly from this mode at small radii) causes — goo 
to rise above unity at sufficiently large radius r (so that this mode there has a 
proper-time frequency that is less than its natural frequency m and hence has a 
concave or exponential radial dependence for larger values of r), one has only the 
asymptotically exponentially decaying behavior of the mode at larger values of r. 

However, the oscillation of the g^ components of the metric (at frequencies that 
are even multiples of m) couples, via the Klein-Gordon equation for the minimally 
coupled massive scalar field, the modes of the scalar field that have frequencies that 
are different odd multiples of m. Thus, for example, the mode with frequency 3m is 
excited and propagates out from the region of the oscillaton where the mode with 
frequency m dominates the stress-energy tensor. 

So far, this is just like the decaying oscillaton described above. If there are 
no incoming waves at frequency 3m, the outgoing waves carry off energy from the 
oscillaton, which decays (and hence is not precisely periodic in time). 

On the other hand, for the precisely periodic but infinite-mass oscillaton, there 
are both incoming and outgoing waves at frequency 3m, and their stress-energy 
tensor (which at sufficiently large r dominates over that of the mode of frequency m, 
since that m-frequency mode is asymptotically exponentially decaying at sufficiently 
large r) eventually causes —goo to rise above 9, so that then the mode with frequency 
3m also has an approximately exponential radial dependence at larger r. With the 
right choice of phase of this mode (i.e., of how big is the part that goes roughly as 
sin(v^8mr)/r at small r), only the decaying exponential is present, and the mode 
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at frequency 3m also goes to zero exponentially rapidly at sufficiently large r. (The 
mode at frequency 3m starts to fall off exponentially in radius at a much larger 
radius than the mode at frequency m does.) That is, the phase is chosen so that at 
each radius where this mode is oscillating radially, there are equal parts of outgoing 
and incoming waves, and then the rising gravitational potential from the stress- 
energy tensor of mostly that mode makes the mode gravitationally bound so that 
the outgoing waves reflect off the potential barrier and become the incoming modes 
at smaller radii. 

Of course, there is also coupling to modes with frequency 5m, and these modes 
must also be given the right phase to have equal parts of outgoing and incoming 
waves for radii at which —goo < 25, so that when one gets so far out that the energy 
density of that mode, mostly, causes — goo to rise above 25 and that mode also 
to develop an exponential behavior, one has only the approximately exponentially 
decaying piece. 

And so it goes for all higher modes with frequencies that are odd multiples of m, 
say (2n + l)m. They must each be chosen so that have equal amounts of outgoing 
and incoming waves for —goo < (2n + 1 ) 2 , where that mode oscillates in the radial 
direction, and then that mode is totally reflected by the gravitational potential for 
-g 00 > (2n + If. 

Since this process must continue indefinitely in order that the oscillaton be pre- 
cisely periodic in time and not be losing energy to outgoing waves of any mode, one 
must have — goo rising indefinitely, so the metric is not asymptotically flat and indeed 
has M(r), the mass interior to radius r, rising indefinitely with r: these are infinite- 
mass oscillatons. (They are also almost certainly unstable to small perturbations, 
but that is another story.) 

These precisely-periodic but infinite-mass oscillatons are of course not physically 
realistic, but they do make interesting theoretical solutions of the coupled Einstein- 
Klein-Gordon equations that one may use for such purposes as defining a precise 
canonical (though rather ad hoc) maximum mass for a finite-mass oscillaton without 
nodes. 

For example, one could take a precisely-periodic oscillaton and then evaluate 
M(t, r) at a time t when the scalar field at that r is passing through zero, and at a 
value of r where the time average of the energy density of the mode with frequency 
m is equal to the time average of the energy density of the mode with frequency 
3m. (If the scalar field passes though zero more than twice in each coordinate time 
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period 2ir/m, then choose one of the times at which the field at that r is zero but 
the magnitude of the time derivative of the field is a minimum among all of these 
times when the field is zero there. If that still does not uniquely specify the time 
within one period of the metric, which is a half-period tt/itl of the scalar field, then 
go to the minimum of the magnitude of the second time derivative of the field among 
those times, etc, until the degeneracy is broken. If the degeneracy never is broken, 
then presumably the mass does not depend on which of those times it is evaluated.) 

The value of M(t, r) at that value of t and of r could then be said to be a 
canonical (though ad hoc) finite precise value for the mass of a decaying finite-mass 
oscillaton without nodes that is in some sense represented by the precisely-periodic 
infinite-mass (at r = oo) oscillaton. 

The idea of the representation is that in the precisely-periodic oscillaton, the 
value of r determined above would be where the mode with frequency m has decayed 
exponentially to a very small value, where its time-averaged energy density has 
dropped to that of the tiny amount of outgoing and incoming waves at frequency 3m. 
For the finite-mass decaying oscillaton, if one does not include the energy density of 
the modes with frequency 3m or higher outside this radius (e.g., if one assumes that 
they are just starting to be emitted and so at that time are negligible outside this 
radius), the energy density of the mode with frequency m is so small at that radius 
(and dropping roughly exponentially with radius) that there is a negligible addition 
to the mass from that energy density in the finite-mass oscillaton. 

That is, it is negligible unless we are wanting to get some absolutely-precise 
number to assign as the mass of the oscillaton, in which case we must go from the 
slightly-poorly-defined initial configuration of a finite-mass decaying oscillaton to 
the precisely defined (but unphysical) infinite-mass oscillaton. 

To get a representation of the maximum mass of a decaying finite-mass oscillaton 
by this sort of precise mathematical definition, one must choose the right precisely- 
periodic oscillaton. There can be an arbitrary number of nodes in the mode with 
frequency m before it decays approximately exponentially in the radial direction, 
and to model the decaying finite-mass oscillaton which has no nodes in its mode 
with frequency m, we want the precisely-periodic oscillaton also to have no nodes 
in that mode. On the other hand, we want the modes at frequencies that are higher 
odd multiples of m, say (2n+ l)m for each positive integer n, to have as low a value 
of energy density possible, which would mean they should each have the largest 
number of nodes possible before —goo rises above (2n + l) 2 [at a rate with respect 
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to radius that is given mainly by the energy density in the outgoing and incoming 
waves at that frequency, once —goo rises above (2n — l) 2 and the mode at the next 
lower frequency, (2n — l)m, is reflected back inward]. 

Smaller numbers of nodes are possible for each of the higher-frequency modes 
if there is an extra magnitude of outgoing and incoming waves at that frequency, 
so there appears to be a whole infinite sequence of integers to be specified for the 
generic precisely-periodic oscillaton, even after specifying, say, the time-averaged 
energy density at the center to give the one continuous parameter of the classical 
oscillaton that is the nonlinear-gravity analogue of the scale of the nearly-Newtonian 
oscillaton. This means that there should be uncountably many oscillatons at the 
same value of the central time-averaged density. 

It is plausible that they form a fractal set of perpetually oscillating spatially in- 
homogeneous but spherically symmetric periodic-in-time solutions of the Einstein- 
Klein-Gordon equation, somewhat analogous to the apparently fractal set of homo- 
geneous but nonperiodic cosmological solutions of the same equations jTTj. 

However, here we want just the simplest example, with the fastest possible falloff 
of the energy density, to use to define a precise value for the maximum mass of 
a finite-mass oscillaton with no nodes for the mode with frequency m and only 
outgoing components for the modes with higher frequencies — these frequencies being 
only approximate in the slightly-non-periodic case in which the finite-mass oscillaton 
is actually slowly decaying with time. 

Once we have chosen the simplest precisely-periodic oscillaton for each possible 
value of the central time-averaged energy density, we still need to choose the one 
that gives the maximum value of the mass defined by the procedure above. That 
will be the unique precisely-periodic oscillaton that we will use to represent in some 
sense the initial state of a finite-mass decaying oscillaton, and it will give us a precise 
(though rather ad hoc) mathematical definition of /x max . 

Now to get a definite initial state [say for giving a precise definition of fj,(t) 
and of its inverse, t(n)] for a finite-mass decaying oscillaton that is represented 
by the unique precisely-periodic oscillaton above, one could take the initial metric 
and scalar field values and time derivatives of the precisely-periodic infinite-mass 
oscillaton at the time t defined above for determining the value of M(t, r) that was 
used to represent the maximum mass of a decaying oscillaton, out to the value of 
r that was also determined by the procedure above (where the time average of the 
energy density of the mode with frequency m had dropped to the same time average 
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of the mode with frequency 3m). One could then truncate the precisely-periodic 
oscillaton at that radius, replacing its initial data on that hypersurface with initial 
data that agreed for smaller values of r but which for larger values of r has zero 
for both the scalar field and its time derivative, and the Schwarzschild values for 
the initial data of the metric. (That is, the solution can be taken to be vacuum 
Schwarzschild on the initial time surface for larger values of r.) 

Next, simply let this truncated oscillaton initial data evolve by solving the 
Einstein-Klein-Gordon equations, to represent the decay of an initially maximal- 
mass oscillaton. 

As energy flows out from the decaying oscillaton in the form of scalar waves 
moving slower than the velocity of light, the ADM mass and also the mass at future 
null infinity stay constant, so we need a different definition of the mass of the 
decaying oscillaton to define a nontrivial time dependence for it. For that, we 
can simply choose M(t,r), with the coordinate value of r kept fixed at the value 
where, in the periodic oscillaton that provided the initial data out to that radius, 
the time-averaged energy density of the mode with frequency m equaled that in the 
mode with frequency 3m. 

Here the coordinate value of r, once determined for the precisely-periodic oscil- 
laton, is assumed to be kept rigid, with no gauge freedom, say by continued use in 
the decaying oscillaton of the gauge choice that goi = and that goo has no periodic 
component. However, since goo has a slow secular evolution for the decaying oscilla- 
ton, I'm not quite sure how to fix the gauge absolutely precisely in this case, since 
I don't see how uniquely to disentangle this slow secular evolution from a periodic 
component. Perhaps if one wanted an absolutely precise mathematical definition, 
which has been the aim of this long discussion, one should use instead as the sphere 
where M(t, r) is to be evaluated as a function of r, a sphere that has a constant cir- 
cumference (the Schwarzschild r-coordinate, after taking r to be the circumference 
divided by 2n, though not the r coordinate that makes goo independent of time for 
the precisely-periodic oscillaton). 

There is also a slight ambiguity in defining t absolutely precisely if goo has a slow 
secular evolution rather than being precisely independent of t as it can be in the 
precisely-periodic oscillaton with a suitable choice for the time and spatial coordi- 
nates. One commonly-used precise way to define it would be to define the constant-t 
hypersurfaces to be orthogonal to the worldlines at constant angular coordinates on 
the spheres of constant circumference, which leads to so-called Schwarzschild coor- 
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dinates. 

One disadvantage of the Schwarzshild (t, r) coordinates is that inside the os- 
cillaton, the worldtubes of constant Schwarzschild-r (circumference divided by 27r) 
oscillate in and out, relative to coordinates in which not only is go% = (as is 
true also for the Schwarzschild coordinates by construction) but also goo is constant 
(for the precisely-periodic oscillaton). Then since the hypersurfaces of constant 
t in Schwarzschild coordinates are by constuction orthogonal to the worldlines of 
constant Schwarzschild r (and of contant angular position on the 2-spheres), these 
hypersurfaces also bend forward and back by a rather considerable amount in the 
interior of the oscillaton, giving periodic effects that depend on this global choice 
of Schwarzschild coordinates and are actually considerably larger than any effect in 
the local geometry. 

However, at the value of r defined above, where the time-average of the energy 
density of the mode with frequency m has dropped down to equal the corresponding 
value for the tiny outgoing and incoming wave modes with frequency 3m, one is 
where the energy density of the oscillaton is so low that the radial oscillations of 
the worldtubes of constant Schwarzschild-r are very small and should not have a 
significant effect. 

One still has to define the t-labelling of this foliation into constant-t hyper- 
surfaces. When the metric is asymptotically flat, as it would be for the decaying 
oscillaton, then one can define it so that t is the proper time along the worldtube of 
infinite circumference. However, with this normalization, the coordinate-time period 
of the approximately-periodic decaying oscillaton would be shifted from the value 
2it/m that it is given by construction (normalization of t) in the precisely-periodic 
oscillaton. 

Another simple but inequivalent choice would be to choose t to be proportional 
to proper time of the central worldline at r = 0. Then one can have the approxi- 
mate period of the approximately-periodic decaying oscillaton be very near 2it/m by 
choosing the constant of proportionality between proper time and coordinate time 
to be the same as it is in the precisely-periodic oscillaton at r = in the gauge for 
that solution in which goo is independent of time. However, one should note that as 
the oscillaton decays and the gravitational potential at the center changes relative 
to that at infinity, the nearly-periodic oscillations of the scalar field will have its 
approximate period of oscillation secularly shifted from the value 27r/m that it has 
initially by this choice of t. 
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A third choice for the labeling of the constant-Schwarzschild-t hypersurfaces 
would be to choose at r = the time coordinate t so that the quantity ip defined by 
Eq. (|2*T|l would be real and positive for all t. This essentially forces the coordinate 
t to be chosen so that, at the center at least, the coordinate period of oscillation of 
the scalar field $ (or of the rescaled dimensionless field 0) is fixed to be 2ii/m. This 
is perhaps the best choice if one wants to count the number of oscillations of the 
scalar field, and to make a precise count even when that number is not an integer. 
Henceforth we shall assume that we have made this choice for the Schwarzschild 
time coordinate t (which will not be the same choice that would make t equal to 
proper time at radial infinity). 

Once a suitably time-evolving sphere is defined to represent the outer surface of 
the decaying oscillaton (with the waves that go outside that radius being considered 
outgoing waves rather than part of the oscillaton), and once one has a precise time 
coordinate t, one can in principle solve numerically for the mass as a function of 
this t. Multiplying the mass M(t) (in time units) by the constant m (the natural 
frequency of the scalar field, in inverse time units) gives fi(t), the dimensionless 
measure of the evolving mass of the oscillaton. 

Assuming that fj,(t) is monotonically decreasing with t (which it conceivably 
need not be within each period, since the waves need not be purely outgoing over 
the entirety of a period, though they should be when averaged over an integral 
number of periods), one can invert this relation to get t(/i) to see the time needed 
for an oscillaton to decay from the maximum-mass configuration to a smaller value 
of \i. In the nearly-Newtonian limit in which \x <C 1, one would expect that this 
time should asymptotically approach that given by Eq. I| 140)1 . not in the sense that 
the difference between the actual t{p) and that given by this formula would go to 
zero, but in the sense that the ratio of the actual to that given by Eq. (|140)l 
would tend to unity as fi tends to zero in this purely classical calculation. 

Of course, it would be interesting to calculate numerically the value of /i max and 
of the ratio of the actual £(//) to that given by Eq. (|140|) . the latter being the 
function 

RQ*) = ^V^tGu). (i4i) 

In particular, it would be interesting to calculate the large dimensionless number 

q = mt(e _ Vmax), (142) 
which is the angle by which the phase of the oscillating scalar field advances (27r times 
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the number of scalar field oscillations) as the oscillaton decays from its maximum 
mass to a mass that is smaller by a factor of e. The numerical calculations done 
above for the nearly-Newtonian limit imply that 



This dimensionless number q, which is somewhat of an analogue of an average 
Q-value of the oscillaton system for masses within a factor of e of the maximum 
mass, is a pure mathematical number determined purely by the Einstein-Klein- 
Gordon equations (with no dependence upon the scale set by the scalar field natural 
frequency m) and by the mathematically-precise (but admittedly rather ad hoc) 
procedure given above. 

It is interesting that q is apparently quite large (at least if the correction factor 
-R(e -1 // max ) is not too many orders of magnitude smaller than unity), a counterex- 
ample to the folklore that numbers defined purely mathematically tend to be within 
a few orders of magnitude of unity. For example, if /i max = 0.607 and if -R(e _1 /imax) 
were unity, then one would get q xs 8.20 x 10 65 . (However, there are other coun- 
terexamples that are even more extreme, such as the unknown prime n that is the 
first positive integer greater than one such that the number of primes less than or 
equal to n, namely 7r(n), is greater than Li(n), the principal value of the integral of 
1/ In a; from x = to x = n, which forms a good asymptotic estimate for 7r(n).) 

It would certainly be worthwhile to calculate numerically the values of /i max and 
of q, as well as the function R(fj) for values of \i up to /i max , but since that cannot be 
done within the nearly-Newtonian calculation reported here, it is beyond the scope 
of this paper and will have to wait for future research. 

It would also be amusing mathematically to get a quantitative description of 
the particular precisely-periodic oscillaton described qualitatively above that I have 
used to represent (by its part within the radius r defined above) the maximum-mass 
decaying oscillaton. In particular, what is the quantitative asymptotic behavior of 
the metric and of the scalar field [each mode of which, of frequency (2n+ l)m, even- 
tually decays approximately exponentially in the radial direction as the integrated 
gravitational effect of this and other modes causes — g 00 to rise above (2n + l) 2 so 
that the mode becomes gravitationally bound]? However, this will also be left to 
future work. 




'max 



(143) 
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10 Quantum Decay of Single-Field Oscillatons 



Besides the classical decay of finite-mass oscillatons to outgoing scalar radiation, 
there are also quantum decay processes that appear to be dominated by the anni- 
hilation of two scalar particles into two gravitons. This rate also goes to zero as \i 
is taken to zero, but only as a power-law in /z, so for sufficiently small /i (depending 
on the ratio of m to the Planck value), it actually dominates over the classical decay 
into scalar radiation. 

The annihilation cross section for two nonrelativistic scalar particles to annihilate 
into two gravitons has been given by DeWitt ^H] to be (where I have inserted the 
factors of the speed of light, c, that he and I usually set equal to unity) 

2nG 2 ml 2n7i 2 G 2 m 2 ,^ At . 
Vnr = = 1 = = , (144 

CrV C'V 

where the first form is in terms of the scalar particle mass in conventional mass 
units, which I have been calling m*, and the second form is in terms of the scalar 
field natural frequency in inverse time units, which I have been calling m, which is 
given by Eq. (J2J) as m = m*c 2 /h. 

Therefore, if we have a number density n of scalar field particles of one species 
(more than one species will be considered in the next Section) , and hence a conven- 
tional mass density 

p = m*n = (hm/c 2 )n, (145) 
the annihilation rate (per time) for one particle passing through is 

R = na NR v = (2ixh 2 G 2 m 2 /c 7 )n. (146) 

Two scalar particles annihilate in each such process, but when one takes the 
square of the number density, there is also a factor of 2 overcounting the number of 
pairs of identical particles, so these two factors of 2 cancel each other and give the 
number rate per time per volume by which scalar particles annihilate as 

' /A Rn = (2nG 2 m 2 Jc 3 )n 2 = (2nh 2 G 2 m 2 /c 7 )n 2 = 2nG 2 p 2 /c 3 . (147) 



dtd{yo\) 

In the nearly-Newtonian limit (which is where the nonrelativistic annihilation 
cross section formula would apply), the mass density with one scalar field present, 
and represented by the dimensionless complex field ip given by Eq. (|21|). is 

m 2 \ib\ 2 

' = -5RT- (148) 
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If one integrates over an oscillaton with one real scalar field (one complex ip), 
the total annihilation rate per time is 



dN 1 



dt 8nc 3 

In particular, for a spherical oscillaton with real i/j, one gets 



rrf\^\ l d i x. (149) 



dN 1 
~dt ~ 2^ 



m ijj r dr. (150) 



Now if we use the fact that by Eq. (j3J), the total mass in time units of the 
oscillaton is M = GM*/c 3 , and if we use the fact that the mass in conventional 
units is the conventional mass m* of each scalar particle times the number N of 
such particles (neglecting the small correction due to the kinetic energy of the scalar 
particles and the gravitational binding energy), M* = m*N = (hm/c 2 )N , then we 
get the rate per time at which the total mass M in time units decays away, the 
dimensionless rate 

dfj, hG dN _ hG f 5 4 3 hG r 5 4 2 



P q = -M quantum = - = R~ m— 77" ~ z — s mipdx = — miprdr. 

d{mt) & dt 87tc 8 J 2c 8 J 

(151) 

For the nearly-Newtonian spherical oscillaton analyzed above, by using Eqs. 
dTo) . (JTfij). P?|l. ((55)1 . (UnnD, and (fTUH|l . one can show that the quantum decay 

rate (per time) in the mass (in time units) is 

P q = "Mquantum « 4^m 2 ^f/i 5 = g^/i 5 = Qt^m 2 /i 5 , (152) 

where 



c 5 



mpi = W— w 1.855 x lO^s" 1 (153) 
V «-G 



is the Planck frequency, 



IhG 



t P} = w 5.391 x 10™ 44 s, (154) 



is the Planck time, the reciprocal of the Planck frequency, and where 
AD 

Q = w 0.008 513 223 934 732 691 876 362 000 723 963 272 14, (155) 

J l oo 

using the numerical results given in Eqs. ()1U2|) and ()104|) . 
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When both the classical decay rate given by Eq. (|135|) and the quantum decay 
rate given by ()152j) are both significant, one has the total decay rate being given by 

It may also be of interest to calculate the expected number of scalarons and 
gravitons emitted during one period of oscillation of the oscillaton scalar field, which 
is a time 2n/m. Since by far the most dominant scalaron emission is at frequency 
3m, the energy of almost every scalaron emitted is 3hm, so the expected number of 
scalarons emitted in one period is 



N _ 2ir/m fdM*c 2 \ _ 2tt /m P A 2 C ^ 

Similarly, the most dominant graviton emission is at frequency m (as two scalarons 
of this frequency annihilate into two gravitons of the same frequency), so the ex- 
pected number of gravitons emitted in one period is 



2n/m (dM^ 
Tim \ dt 



quantum 

,5 



w 0.053 490 163 543 442 036 463 830 655 145 560 622 6 /i 5 . (158) 

It is interesting that although it is the classical emission power (into scalar waves 
or scalarons) P c that depends only on fi, with the quantum emission power (into 
gravitons) depending also on m/rapi, for the expected number of particles emitted 
in one oscillaton period, it is the quantum emission into gravitons that depends only 
on fi (with an expected number of gravitons per period never larger than unity, and 
in fact never larger than roughly 0.0044 if the maximum value for /x is 0.607 and if 
the formula above indeed applies to this large a value of //), whereas the emission 
into scalarons (classical) gives a number depending also on m/mpi (and which can 
be larger than unity for sufficiently small mass m). 

The expected number of scalarons emitted per oscillaton period is unity, N s = 1, 

for 



e 2 ^ 

mpi yU 

2 820.165 789 522 802 746 688 836 089 134 773 52 i^bb 897 598 sso cm 546 85 9 reo 732 713 1 84 8 
~ e c 

(159) 
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If this formula were to apply for \i = 0.607, which is one numerical result for the 
maximum mass parameter of an oscillaton [8 , then this would occur for 

m* « 3.63 x 10" n m, P1 « 4.43 x 10 8 GeV/c 2 . (160) 

Then there would be some allowed oscillaton mass value at which the oscillaton 
would emit an expected number of one scalaron per oscillation period, for any 
scalaron mass less than some scalaron mass value that is given crudely by Eq. ()160|) 
that uses Eq. ()135j) outside its jj, <C 1 domain of validity. 

If we define for use here (not the same x as used previously for the rescaled radial 
variable) 

a 39.433 795197160163 093 719 521465 426 369 6 .,„.,, 
X = — RS (161) 

fl fl 

and 

C 

7 = — - « 0.003 008 268 339 955 585 528 843 272 770 555 262 66, (162) 
a'Q 

then the ratio of the classical mass-loss rate to the quantum mass-loss rate is, for 
M< 1, 

# = 7^*V, (163) 
_r g m 2 



which is unity for 



rn 



2 

PI 



x-71na; = lnl7-^l, (164) 



or 



x = i n 7 !> +71nx = ln +71n In +71nx 

y m l ) \ m 2 J \ \ m 2 J J 

= ln ( 7 ^) +71n ( hl ( 7 ^) +71n ( ln ( 7 ^) +7111I )), (1M) 

to give the first three steps in an iterative procedure for solving for the value of 
x = a/ fi that gives P c = P q . 

The decay time for an oscillaton to have its dimensionless decay parameter decay 
from /ii to a smaller value of \x is given by 

m(t- tl )= r-p^p- (166) 

Define the m-dependent constant r [for the use in this immediate section; not 
the same r as used elsewhere for the radial variable or with yet another meaning 
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as the dependent variable defined in Eq. ()123j) ] as the largest real solution to the 
equation 

r | 

S = ( 167 ) 



(lnr) 3 c^Q^m 2 
where the m-dependent constant s here is also not the dependent variable defined 
in Eq. (|124|) . For brevity also define here q [now neither the dependent variable 
defined in Eq. ()122j) nor the large dimensionless constant angle defined in Eq. (|142|) ] 
to be the function of the classically dimensionless mass /i that is 

(168) 

Then when Ins = — In (^Qt^m 2 ) ^> 1 (which implies that lnr 3> 1), and when 
Eq. ()156j) also holds, one can show that the integral on the right hand side of Eq. 
(|166|) can be approximated by the following explicit function of r, q, and q\ (the 
value of q when \i has its initial value 

m(t - h) « {[In (r + g)] 4 - [In (r + gi )] 4 }. (169) 

We can use Eq. (I167j) to rewrite this equation in the form 

Aa 4 Qt 2 Pl m 3 (t - h) « [In (r + g)] 4 - [In (r + gi )] 4 . (170) 

This gives an approximate formula for the decay time, t — 1\, for an oscillaton 
in terms of its mass M and the scalar field mass m (which then determine the 
dimensionless mass parameter /i), assuming that one knows the initial value \i\ of 
the dimensionless mass parameter. 

For other purposes, one might know the decay time and from it want to get a 
relationship between the oscillaton mass M and the dimensionless mass parameter \i. 
For this purpose, it is helpful to define from the decay time the large dimensionless 
parameter [not the same a as the //-dependent constant defined in Eq. ([118)1 ] 

/ 16 \ 1/3 /t-tA 2 / 3 , , 

° s [*q) (V) >>L (m) 

It is also helpful to define a new dimensionless mass parameter 

M 

V= [A^QtUt-t^ (172) 

which depends on the oscillaton mass M and on the decay time t — 1\ but not on 
the mass m of the scalar field as the other dimensionless mass parameter \x = Mm 
does. 
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Then, for example, Eq. (jl70j) becomes 



fi « u{[\n (r + g)] 4 - [In (r + <?i)] 4 } 1/3 , (173) 
and Eq. ()167|) becomes 



av 2 



(lnr) 3 /i 2 
For s ^> 1, this solution to this equation is very roughly 



(174) 



r . s(lns) 3 = _^_j . (175) 

Now we can insert r from Eq. ()175|) and the definition of q from Eq. ()168j) into 
Eq. ()173j) to get the following rough explicit algebraic relation between /i, fix, and 
for fixed t — ti, and hence for fixed a defined by Eq. (|171j) : 



In 



fi 2 \ n 2 J aC 



In 



fi 2 \ fi 2 I aC 



(176) 

However, I should note that although this relation is a good approximation to Eq. 
()173j) when /z is very small, it can be off by about 20 — 30% when /i is large. 

For a/fi\ < a/fi <C In a, so that the classical decay dominates, this relationship 
can be solved explicitly for v to give 

v w (177) 

4(/i 4 e a /^ -/ife Q /^)' 

In this limit, this Eq. (jl77|) is actually a better approximation to Eq. ()173|) than is 
Eq. ()176|) . At the other extreme, for a/fi ^> In a, so that the quantum annihilation 
dominates during most of the decay, this relationship can also be solved explicitly 
for v to give a fairly accurate approximation to Eq. f)173j) . namely 

(4 4 \ -I/ 3 

l_^_ 4 ^ln^- . (178) 
Hi a aC J 

In between these two extremes, i.e., for a/fi ~ In a, I don't see how to give any 
simple expression that would solve explicitly for either [i or v in terms of the other 
(say for fixed a and /ii), though of course one could solve Eqs. ()173|) and ()174|) 
numerically for either /i or v in terms of a specific value of the other, for fixed decay 
time t — t\ and hence for fixed a given by Eq. (j!71|) . 
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Let us put in some numbers for these quantities. If we take fix to be the maximum 
mass of an oscillaton given by [Hj, namely \i\ = 0.607, then we get 

u 4 

qi = r±-e a/tM1 « 1.48 x 10 19 . (179) 
aC 

If we take the oscillaton decay time t — t\ to be the present age of the universe, 
about 13.7 billion years or 8.02 x 10 6 %i in Planck units, then 

1 6 /f _ iA 2/3 

« 3.68 x 10 39 , (180) 



or In a f« 91.10. 

In the previous paragraph it was noted that Eq. ()177|) applies for a/fi <C In a 
and that Eq. ()178|) applies for aj\i ^> In a, so these inequalities are saturated and 
one has a/fi = In a at /i ~ 0.433. Therefore, since \x cannot be much larger than 
this value, one never really has the validity of the inequality a/fi <C In a, but in 
actuality Eq. (|177|) is a good approximation to Eq. (|173j) for e a ^ 1 <C a, and indeed 
at \i = fx u one has e a/ ^/a « 4.44 x 10 -12 < I. 

However, one must still remember that Eq. ()173|) . or equivalently Eq. ()169|) or 
Eq. ()170|) . are valid approximations only to the extent that Eqs. ()135|) and (|152jl 
are valid for the classical and quantum decays rates P c and P q respectively. In this 
paper these formulas were derived under the assumption that \i <C 1, so they are 
not likely to be accurate for the small \i values where e a ^ a, particularly for 
P c with its very strong dependence on fi. Hence the formulas used here for when 
the classical decay appears to dominate should be taken with a big grain of salt, 
as merely provisional formulas that might give some rough qualitative indication of 
the true quantitative behavior. 

Nevertheless, to get some idea of this rough qualitative behavior, let us for now 
assume that Eqs. (j!35|) and ()152j) are valid for all fi less than its maximum value, 
at, say, \i = \i\ = 0.607, and put in some various possible values for m, /i, and/or 
M. 

For example, taking the example used by Seidel and Suen |lj in which the scalar 
field mass is typical of that of an axion, m* = 10 _5 eV/c 2 or m m 1.519 x 10 10 s _1 ~ 
8.19 x 10~ 34 m P i, then ln^p^/m 2 ) w 146.564, so one finds that P c = P q at x ss 
183.031, or n w 0.2154 (if the formulas above really apply to this large a value of 
/i), which corresponds to an oscillaton mass in time units of M ~ 1.418 x 10 _11 s or 
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an oscillaton mass in conventional units of M* « 5.726 x 10 27 ^ = 0.9585M© (about 
96% of the mass of the earth). 

For this oscillaton, assuming that the equations above did apply for \i pa 0.2154 
even though this is not much smaller than unity, one gets that the total power 
emitted would be P = P c + P q ~ 5.3 x 10~ 72 , and the logarithmic rate of decrease of 
the mass would be —d\nM/dt pa 3.7 x 10 s -1 . The oscillaton would contain about 
iV = M*/m* = (mpi/m) 2 /i pa 3.2 x 10 65 scalarons, and so in each period of oscillation 
of the oscillaton, there would be about (2ir/3)(mpi/m) 2 P c ~ 8.3 x 10~ 6 scalarons 
emitted (each of energy roughly 3m*c 2 ) and about (2n)(m-p\/m) 2 P q ^ 2.5 x 10 -5 
gravitons emitted (each of energy roughly m*c 2 ). That is, one would need to wait 
on average about 120 000 periods of oscillation between the emission of successive 
scalarons, and about 80 000 periods of oscillation between the emission of successive 
pairs of gravitons (since they come out predominantly in pairs, with the pair having 
two-thirds the energy of a typical scalaron that is emitted). Therefore, although 
this oscillaton is not absolutely stable, for astronomical purposes it is very nearly 
stable. 

If this oscillaton with /x pa 0.2154 had actually decayed from /i = \i\ — 0.607, 
that would have taken a time t — t\ ~ 10 51 yr, again using Eqs. (|135|) and ()152j) 
outside their true range of validity just to give a qualitative answer. 

On the other hand, if an oscillaton with this value of the scalar field mass, 
m* = 10 _5 eV/c 2 , were to have started at \i\ = 0.607 at the beginning of the universe 
and had decayed up until its present age, it would now have \i ~ 0.459 (again taking 
the classical decay rate formula outside its range of validity, /i < 1). But even 
though this calculated value for \i today is not likely to be actually correct, it is 
interesting that it is significantly below the initial value. Thus if this result is at 
least qualitatively correct, an initially maximum- mass oscillaton with this value of 
the scalar field mass would have decayed by a significant amount during a time 
comparable to the age of the universe. 

To take a more extreme example, if we imagine that there is a scalar field 
(quintessence?) of natural frequency m that has the value of the current Hub- 
ble expansion rate, Hq pa 2.3 x 10 _18 s _1 ~ 1.24 x 10~ 61 mpi, which corresponds to 
m,c 2 pa 1.51 x 10~ 33 eV, then In inlJm 2 ) pa 274.678, so then P c = P q pa 4 x 10~ 129 
at x pa 314.945, or jj, pa 0.1252, giving an oscillaton mass in time units of M ~ 
5.4 x 10 16 s ~ 1.7 billion years, or an oscillaton mass in conventional units of 
M* pa 2.2 x 10 55 (? = 1.1 x 1O 22 M (of the same order of magnitude as that of 
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all the observable galaxies in the universe, by the same coincidence that this mass 
is roughly within a factor of 10 or so of what is needed to close the universe). For 
this example, one would need to wait on average about 1 800 000 periods of oscil- 
lation between the emission of successive scalarons, and about 1 200 000 periods of 
oscillation between the emission of successive pairs of gravitons. 

If we took an oscillaton with a scalar field of this m = Ho, started it with 
[L\ = 0.607, and applied the formulas above, we would find that during the age of 
the universe, it would have had fi decay only by about 6.3 x 10~ 22 , an insignificant 
reduction in its value. Besides the usual caveat about the inapplicability of Eq. 
()139j) to this large \x value — note that Eq. (|140J) does not apply here, since in this 
case the two terms in Eq. (j!39|) are very nearly equal — there is also the error from 
the fact that the age of the universe corresponds to only about 0.99 of a radian of 
the phase of the oscillation of an oscillaton with a scalar field mass equalling the 
current value of the Hubble constant, whereas the formulas above apply only for 
time periods containing many oscillations. 

For values of x smaller than the solution of Eq. (j!64j) or ()165j) . so that fi is larger 
than the corresponding critical /i value for that m/mpi, then the classical decay rate 
dominates (i.e., for large oscillaton masses). On the other hand, for values of x larger 
than the solution of Eq. (j!64J) or (jl65J) . so that [i is smaller than the corresponding 
critical /i value for that m/mp\, then the quantum decay rate dominates (i.e., for 
small oscillaton masses). 

It is interesting that even with perhaps about the smallest conceivable value 
of the scalar field mass in the present universe, that of the Hubble constant, the 
quantum emission dominates over the classical emission when fi is only as small 
as about 1/8 (and of course for all smaller values of fi). That is, for almost any 
conceivable oscillaton in the present universe, if the dimensionless mass parameter 
/i is smaller than roughly 1/8, the classical emission of scalar waves would be even 
less than the tiny quantum emission of gravitons from the annihilation of pairs of 
scalar particles in the oscillaton. This illustrates how rapidly the classical emission 
drops as \i is made small. 

When the mass of the scalar field is much smaller than the Planck mass, as it 
must be for one to have nearly-Newtonian oscillatons containing a large number of 
scalar particles, as is implicitly assumed in the analysis above, then by Eq. ()164|) 
or (|165jl one finds that x is large in comparison with unity for P c = P q . However, 
because a is itself rather large, it is not necessarily the case that this value of x 
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corresponds to a small value of fi or even a value of fi less than /i max . 

For example, if Eq. ()163j) were valid not just for very small fi but also for 
values of \i up to /i max , and if we take /i max = 0.607, then we would find that 
P c = P q at fi = /i max = 0.607 for m ps 9.47 x Kr 10 m P1 or m* ps 9.47 x KT 10 m* P i ps 
1.16 x 10 10 GeV/c 2 . Assuming that this is correct, then for larger values of the scalar 
field mass, the quantum decay would dominate (P q > P c ) for all values of \i up to 
and including /i max . On the other hand, for smaller values of the scalar field mass, 
which is more realistic if the scalar field is, say, an axion, then there is always a mass 
range (say \i c < fi < ^ max with fi c being given roughly by a/x with x the solution of 
Eq. ()164j) or ()165j) when it gives a small fi c ), where the classical emission dominates 
for the mass loss rate, though for fi < fi c < // max the quantum mass-loss decay rate 
would dominate. 

It may be of interest to estimate the present upper bound on fi for oscillatons with 
scalar field masses other than the two examples above, assuming that the oscillatons 
formed in the early universe and have been decaying for a time comparable to the 
age of the universe. The maximum present value of fi that they would have would 
be what they would have if they started with the maximum allowed initial value of 
fi, which here I shall take to be fi\ = 0.607, as above. I shall also assume that the 
scalar field mass-energy is m^c 2 <C 10 10 GeV, so that the decay from fi = fi\ within 
the present lifetime of the universe would be in the regime where the classical decay 
dominates (see below for more details on this) and where I shall assume that Eq. 
(HH holds. 

Then if I use x = a/fi defined by Eq. (jlfilj) . one gets (with a new use for y) 



x 

where for this Section I shall define 
C 

V = -n[m(t - ti) + qi] 



e 



;isi) 



m Jf c 2 \ ( t — U 



w 4.067 595 xl0 34 — ^ ^7— + 9.188 619 x 10 20 (182) 

V eV j \ L37 x 10 10 yr/ 

(not the same y that denoted the rescaled imaginary radial coordinate in Section 8). 

An explicit approximation that solves Eq. (|181|) for x in terms of y to at least 
8-digit accuracy is 

x ps In y + 4 In (lny + 4 In (lny + 4 In (In 7/ + 4 In (In 7/ + 41nxi)))), (183) 
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where 

ol 

Xl = _ ~ 64.965, (184) 
ilnxx = 4rn(o;//x 1 ) « 16.695. (185) 

Then we get 

39.433 795 , , 

n ~ (186) 

ln?/ + 41n(ln?/ + 41n(ln?/ + 41n(lny + 41n(ln?/ + 16.695))))' 

However, Eq. f)181|) itself is not that accurate, since it was derived on the as- 
sumption of the accuracy of Eq. ()139|) . which is in doubt, since that equation was 
derived for /i < 1, whereas here for m*c 2 <C 10 10 GeV, one gets fi in the range 
roughly between 0.3 and 0.6, which is not much less than unity. 

Once we have an estimate for fi (no doubt rather crude, since it does not give 
<C 1 where it would be valid), or for the maximum value of function 
of the scalar field mass, we can easily get the oscillaton mass M = fi/m in time 
units or M* = hcfi/(Gm*) in conventional mass units. In terms of the solar mass 
M « 1.989 x 10 33 g « 0.9137 x 10 38 m^ P1 , one can use Eq. (gj) to write 

M* « 1.336 337 63 x 1O~ 1O M (^-^) \i. (187) 

Combining this with Eq. (J186)) then gives 

5.269 686 43 x 1Q- 9 M [(1 eV)/(m*c 2 )] 



In y + 4 In (In ?/ + 4 In (In y + 4 In (my + 4 In (In ?/ + 16.695)))) 

where y is given by Eq. (|182j) . This would be the estimated value for the oscillaton 
mass if it started at fi\ = 0.607 and would be an upper limit for the mass if 0.607 
were the maximum value of /i at which it could have started. Again, this formula is 
applicable for m^c 2 -C 10 19 eV (where the classical decay dominates for an oscillaton 
starting with Hi = 0.607 and decaying for up to 13.7 billion years), and it also 
assumes the dubious correctness of Eq. ()135|) for the resulting fairly large values of 
/i as given by Eq. (|186j) . 

For example, if we let fj,(m*c 2 ) be the value of /i that an oscillaton of of scalaron 
mass-energy m*c 2 would decay to, from \i\ = 0.607, in a time of 13.7 billion years, 
then we had shown above that /i(10 _5 eV) ~ 0.459. We can also readily calculate 
the following values of /j,(m*c 2 ) for other values of m*c 2 : 

/i(10~ 35 eV) w fix - 4.41 x 10" 24 , 



59 



/i(l(T 30 eV) 


« /ii - 4.41 x 1(T 19 


/i(10- 25 eV) 


« /i! - 4.41 x 1(T 14 


/i(10- 20 eV) 


ss /ii -4.41 x 1(T 9 , 


/i(10" 15 eV) 


« /ii -4.41 x 1(T 4 , 


/i(10- 10 eV) 


« 0.534, 


/x(l(T 5 eV) 


« 0.459, 


/i(l eV) 


« 0.402, 


//(10 5 eV) 


« 0.358, 


/i(10 10 eV) 


« 0.323, 


/i(10 15 eV) 


« 0.295. 



(189) 

Similarly, we can calculate M*(m*c 2 ), the value of M* that an oscillaton of 
scalaron mass-energy m*c 2 would decay to, from ji\ = 0.607 and 



M„! = ^ » 8.111 x lO- n M i^, (190) 



in a time of 13.7 billion years, as having the following values: 

/i(10" 35 eV) « 8.111 x 1O 24 M , 

//(10- 30 eV) « 8.111 x 1O 19 M , 

/i(10" 25 eV) « 8.111 x 1O 14 M , 

/i(10~ 20 eV) « 8.111 x 1O 9 M , 

/x(10" 15 eV) « 81O55M , 

/i(10~ 10 eV) « O.7133M , 

//(10 _5 eV) « 6.128 x 1O" 6 M , 

/i(l eV) « 5.375 x 1O _11 M , 

/i(10 5 eV) « 4.790 x 1O~ 16 M , 

/i(10 10 eV) « 4.322 x 1O~ 21 M , 

/i(10 15 eV) « 3.938 x 1O~ 26 M . (191) 

These are either the estimates for the masses, if the oscillatons started 13.7 billion 
years ago with \i = 0.607, or are estimates for the upper bounds of the oscillaton 
masses, if 0.607 is merely an upper bound on the initial value of /i. 
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11 Quantum Decay When the Particle Number 
Gets Small 



When fi <C 1 and when Eq. ()163|) gives P c /P q <C 1, then the quantum mass-loss 
rate dominates and is given to good accuracy by Eq. (jl52J) . but only so long as the 
number of scalar particles in the oscillaton, 

JV«^ = ^=-^L = I^« (192) 

is large in comparison with unity. 

That is, the quantum mass loss rate dominates and is given to good accuracy by 
Eq. (HS2D if 

m2 </i</i~ a - ~ 19 - 716 897 598 580 082 < 1 r 193 \ 

ra P1 x ln(7mp 1 /m 2 ) 61.769 + In (eV/m*c 2 ) 

The right hand side of this requirement is actually a bit stronger than what is 
needed, which is that both /i < \i c and that // < 1, but it is not really necessary 
that /i c 1. We may note that for this inequality to have any range of validity for 
/i, we need m <C mpi, which we have been assuming thoughout this paper and shall 
continue to assume. 

When an oscillaton is decaying, it will eventually get down to having a small 
number N of scalar particles, and Eq. ()152j) will cease to be accurate. In principle 
one could solve the iV-body Schrodinger equation with Newtonian attractive poten- 
tials between the N scalar particles for the ground-state wavefunction (ignoring for 
the moment the annihilations into gravitons) and then calculate the overlap between 
two particles to get the two-particle annihilation rate into two gravitons. However, I 
did not do this calculation for N > 2 and am not familiar with the literature where 
it might have been done. 

Just as the annihilation rate for large N goes as the 5th power of \x and of N, 
one would also expect that the annihilation rate for a small number iV would also 
decrease rapidly as N is reduced, reaching a minimum for N — 2 (if one can reach 
this number, though if one has an odd number when N is somewhat larger, and if the 
2-body annihilations dominate so that the scalar particles predominately annihilate 
in pairs, then one would most likely end up with a 3-particle state before the final 
decay to two gravitons and one free scalar particle). 
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So if the oscillaton decays down to two scalar particles before annihilating com- 
pletely, the decay of the final two particles is likely to take more time than the entire 
decay down to that point. 

The final annihilation rate is easily calculable from using the ground-state solu- 
tion of the two-particle Schrodinger equation. One readily gets that the probability 
density for one of the particles to be at the location of the other, the quantity that 
takes the role of the number density n in Eq. (j!46|) . is 



n = A 6-3' ( 194 ) 



continuing to use m with units of inverse time. Then by Eq. (j!46|) . the annihilation 
rate per time (for the two particles to annihilate) is 

R= T?T . (195) 

If one uses the first Eq. (J151)) to convert this to the quantum expectation value 
of a mass-loss rate and uses the fact that for this 2-particle state, fi = 2(m/mp{) 2 , 
one gets 

dM d\\i 



Pq — -^quantum 



dt d(mt) 



m \ 12 h 6 G 6 12 1 m 2 

-m 12 = —^-pL 



mpi / c 30 32 m P1 

« 3.670 759 777 914 995 479 155 135 882 831 790 22 gt P1 mV- (196) 

That is, the actual rate at which the 2-particle state annihilates is a factor of about 
3.67 times what one would get by blindly extrapolating down to N = 2 particles 
the rate given by Eq. ()152|) . which actually applies only for very many particles, 
iV ^> 1 (as well as ji <C 1). Eq. (|196|) contains the largest positive power of ft (6), 
the largest positive power of G (6), and the largest negative power of the speed of 
c (30), that I can ever recall seeing in a formula, though I am not used to using 
formulas in which I have not just set h — G — c — 1. 

If instead we take the reciprocal of the annihilation rate R given by Eq. (J195)) 
as the expectation value of the decay time from the 2-scalar-particle state to the 
2-graviton state, and also as an estimate for the total decay time for an oscillaton 
(since it presumably dominates over the time to get down to 2 particles, assuming 
that the number of particles is even when one gets close enough to the 2-particle 
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state that one can ignore the probability that an odd number of scalar particles will 
annihilate), then we get a total decay time of 

2r4° 2c 25 2h 6 c 3 /m.pA 11 

t decay sa l/R = — — = = - 2% — - . (197) 



m 11 K b G*m xl G^m\ l ri V m* 

For example, if we take a typical axion mass, = 10~ 5 eV/c 2 , then Eq. f|19Cj|) 
gives P q w 9.11 X 10 -398 and Eq. flTTJTI) gives t d eca y = 1-88 x 10 346 years. To take the 
more extreme example in which m has the value of the current Hubble expansion 
rate, H « 2.3 x lO -1 ^ -1 « 1.24 x 10 _fll m P i, then Eq. ffM gives P g « 1.3 x 1(T 731 
and Eq. (fTTT7jl gives t de ca y ~ 1.3 x 10 680 years. 

Thus the complete quantum decay of an oscillaton can take a very long time 
and probably would not be a suitable subject for an experimental Ph.D. thesis. On 
the other hand, the slowness of both the classical and quantum decay of oscilla- 
tons of light scalar fields shows that if they form in the universe, they can last for 
astronomically long times. 



12 Quantum Decay of Multiple-Field Oscillatons 

The analysis above was for oscillatons having just one real massive scalar field, min- 
imally coupled to gravity. Although in this paper we shall not go beyond minimal 
coupling or consider other scalar- field self couplings (other than mass terms), we 
started with a general discussion of an arbitrary number of minimally coupled mas- 
sive scalar fields and their classical decay rates, so it would be of interest to say also 
how the quantum decay goes when there are more than one scalar field. 

When the scalar fields all have different masses, then the separate decay processes 
are incoherent, so the rates for each just add, with the rate for each (in the nearly- 
Newtonian limit) going as the spatial integral of the square of the mass density for 
that scalar field, with the coefficient as given above. 

This is at least so if we average over times long in comparison with the reciprocals 
of the differences of the scalar field masses in frequency units, which will hereby be 
assumed — if any mass differences are short in comparison with the reciprocal of the 
decay time of interest, then for that time we may consider these scalar fields as 
having the same mass. Intermediate cases in which the decay times of interest are 
comparable to the reciprocal of any mass differences will not be discussed here. 

Therefore, we may consider separately all of the fields at one mass (or one range 
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of masses if the range is much less than the reciprocal of the decay time being 
considered) . 

The simplest case is that in which there are two equal-mass scalar fields that are 
oscillating at 90° out of phase. This is equivalent to one complex scalar field that 
has a global U(l) symmetry and hence a conserved particle number that presumably 
cannot decay away, at least by perturbative quantum effects such as what DeWitt 
[T5] used to calculate the annihilation of scalar particles into gravitons. 

Presumably there are nonperturbative gravitational effects in which a nonzero 
particle number, though conserved by the global U(l) symmetry perturbatively, 
forms or tunnels into a real or virtual black hole that then decays into a different 
particle number (e.g., zero). Thus at some level the global U(l) invariance is surely 
broken by gravity. In ^Hj we used a model of gravitational foam to estimate that this 
rate would be disastrously high for point scalar particles, suggesting that perhaps 
no such particles could exist in our universe. If so, this would of course rule out 
the whole idea of oscillatons (unless they were made of composite scalars that are 
not pointlike down to near the Planck scale). But since our ideas were admittedly 
rather speculative, here I shall assume that the nonperturbative effects violating 
global U(l) invariance are suppressed to give rates much smaller than the particle- 
antiparticle annihilation into gravitons that is allowed by the perturbative analysis 
that preserves the global U(l) invariance. 

If this is indeed so, when there are two equal- mass scalar fields that are oscil- 
lating at 90° out of phase, effectively the decay of each individual scalar field must 
destructively interfere so that the total decay rate is zero. 

At first this sounds impossible, since if one has a state with N\ scalar particles 
of the first field and N 2 of the second, then the final state in which two particles of 
the first field annihilate into gravitons would have N% — 2 particles of the first field 
and N2 of the second, which would be orthogonal to the final state in which instead 
two particles of the second field annihilate into gravitons, leaving Ni particles of 
the first field and N 2 — 2 particles of the second field. Therefore, how could there 
possibly be any destructive interference to prevent the particles from annihilating? 

However, this objection can be circumvented if the quantum state of the oscilla- 
ton does not have a definite number of particles of both kinds. (Indeed, that would 
have to be the case if they are oscillating 90° out of phase, since phase is in some 
sense a conjugate variable to particle number. Note that the total particle num- 
ber could be precise, so that the total phase is undefined, so long as the individual 
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particle numbers are sufficiently indefinite that the relative phase between the two 
fields is well defined.) For example, the quantum state for the particles could be a 
coherent state that is an eigenstate of the annihilation operators for the two kinds 
of particles. 

Then in the case that two particles of the first real scalar field decay, although the 
expectation value of the number of particles of that field would have been reduced 
by two, the final state would not need to be orthogonal to the state that would result 
if instead two particles of the second real scalar field were to annihilate. Therefore, 
the two decay processes can interfere. When the two scalar fields oscillate 90° out of 
phase, their combination is equivalent to a single complex scalar field with a U(l) 
symmetry that prevents perturbative quantum decay into gravitons. 

The gravitational signal of this U(l) symmetry would be that the stress-energy 
tensor would have no oscillations, so one would be back to the case of a boson star 
that seems to be completely stable (except presumably to nonperturbative tunneling 
processes in which some or all of the particles tunnel into a black hole, or a virtual 
black hole, that would either transcend or violate what would otherwise be the 
conservation of the global U(l) charge [19]). 

Now consider the case in which there are an arbitrary number, say n, scalar 
fields at some mass mj. We shall continue to assume the nearly-Newtonian limit, in 
which the dimensionless rescaled real massive scalar fields <f>u, defined by Eq. (jSjl . 
have the form given by Eq. (|17jl in terms of the complex dimensionless scalar fields 
tpH that are very slowly varying spatially, and even more slowly varying temporally, 
if at all, on the length scale c/mj and on the time scale 1/mj. Then the real scalar 
fields (fiij are essentially oscillating nearly periodically with frequency mi. 

When two scalar particles annihilate into (predominantly) two gravitons, the 
graviton wavelengths are roughly c/mj, which is a much shorter length scale than 
the length scale of the variation of the fields <pu, so each region of size of the order 
of c/mj annihilates essentially independently. (That is, the quantum states of the 
outgoing gravitons are essentially orthogonal for the annihilation in the separate 
regions, if one uses graviton wavepackets that have sizes more nearly comparable to 
their wavelengths than to the much bigger size of the oscillaton.) Thus we can just 
add up the annihilation rates in each region, effectively getting an integral over the 
oscillaton of the annihilation rate in each region. 

Since we shall only be interested in the annihilation over many oscillations of 
the oscillaton, we shall only consider the annihilation rate averaged over many such 
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periods. 

In each region of size somewhat bigger than c/mj where we are calculating the 
average annihilation rate, each of the real scalar fields of mass mj is oscillating 
essentially with constant amplitude and period, staying in phase with each other 
scalar field over a time long compared with the oscillation period that is very nearly 
27r/m/. In this region, we can perform an 0(n) transformation of the n scalar fields 
so that all but two of the fields are transformed to zero for the time of interest 
(to the accuracy of the nearly-Newtonian approximation), and the two that remain 
nonzero are oscillating 90° out of phase, say 

0i = E OiJ^u ~ 2c i cos M* - *o)], (198) 
j 

02 = E °2J0/j « -2c 2 sin [mj(t - *„)], (199) 
J 

and 

& = £CWjj«0 (200) 
j 

for i > 2, with real positive amplitudes C\ and c 2 that are nearly constant over 
the spacetime region where the time-averaged annihilation rate is being calculated. 
Without loss of generality we can choose the 0(n) transformation to give c\ > c 2 . 

I am not bothering to include the subscript I, which tells what the mass m I is, 
on what I am calling 0j. Indeed, the fact that I am giving 0j only a single subscript 
is used here to distinguish it from the 0(n)-related scalar fields 0/j without having 
to put primes on 0j as I would have if it had the same number of indices as 0/j. 

We can also define the 0(n) -transformed complex scalar fields 

^i=E°i^/J~c ie ^ t0 , (201) 
j 

^2 = E « -ic 2 e miU \ (202) 

j 

and 

^ = E°^/j~ (203) 
j 

for % > 2. 

Now for the quantum analysis, we can replace the two real scalar fields with the 
one complex scalar field 

$ = ^(0 1 + ^ 2 )^^[(^ 1+ ^ 2 )e^+(^ 1+ # 2 )e^] 

w 4=[(ci + c 2 )e- im/( *-* o) + (ci - c 2 )e imi{t - to) ] = <S> + e- imit + $_e imjt (204) 
v2 
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(with $ not to be confused with the original real scalar fields $jj or the slowly- 
varying complex dimensionless scalar fields ipu that were used to represent each 
dimensionless real scalar field <fiu] those all had subscripts that will not appear on 
the single rapidly-varying complex $ that combines the two rapidly-varying real 
scalar fields 4>\ and ^2)- 

Because of the [/(l)-invariance of the field equations and stress-energy tensor of 
the two nonzero real scalar fields, 0i and 02, which are represented by this complex 
scalar field $, there is a conserved global U(l) charge. The part of the classical field 
$ that has the phase factor e -* m /(*-*o) anc i fa e coefficient <3>+ can be said to represent 
particles with positive global U(l) charge and with particle number density 



n I+ 



mjc 
AnGh 



mic 
AnGh 



ci + c 2 



2V2 



mjc 
AnGh 



V2 



(205) 



and the part of the classical field $ that has the phase factor e +mi ^~ to ^ and the 
coefficient $ + can be said to represent antiparticles with negative global U(l) charge 
and with antiparticle number density 



ni- 



mjc 
AtxGTi 



l$-l 



AtxGTi 



Ci - c 2 



2^/2 



mjc 
AnGfi 



V2 



(206) 



In the classical limit that we are assuming, we can express these number densities 
in terms of the mass density and mean-squared pressure of the fields of mass mj in 
the following way: 

By using Eq. (J26|) for Goo ~ 8irGp or Eq. pTjl directly for the mass density p 
(which is nearly constant in time), and splitting it up into the contributions from 
the fields of the different masses mi, one gets 



Pi 



where 



m 



u\ 



4ttG 



c +c 



mi% 



-m 



Hi = n I+ + 



mjc 



mjcr 
AnGh 



1 j i 



(207) 
(208) 



AnGh y 1 z 

is the total number density of all the scalar fields of mass mi, or, equivalently, of 
both the particles and the antiparticles of the single complex scalar field $ that 
classically represents all of the real fields. (I.e., we are ignoring vacuum fluctuations 
in the transformed scalar fields <f>i with i > 2 that are classically zero.) 
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Similarly, by using Eq. (|25jl for Gy ~ (87rG/c 2 )P5ij for the oscillating nearly- 
isotropic pressure P, and also splitting it up into the contributions from the fields 
of the different masses mj, one gets 

Pi^-IT^ E(^e- 2im/t + ^e 2 ^). (209) 

o7T(_r j 

In this case, if one takes the time-average of the square of the total pressure of all 
the n scalar fields of mass mi (or equivalently of the single complex scalar field $), 
one gets 

* TVi^J ( c ? - c 2f ~ 2 ^V+n/- (210) 

Because P/ is oscillating sinusoidally, its maximum value, say P/ max (as a function 
of time at each spatial location) is (2P 2 ) 1 / 2 . 

Then from Eqs. (|2()7jl . (|2()8jl . and (|21()|1 . one can solve for the number densities 
of both the particles and the antiparticles of the complex scalar field <3> of mass mf. 

where (2P 2 ) is given in terms of the ipi/s by Eq. (|210j) . 

Another way to express this relationship, using P/ max = (2P 2 ) 1 / 2 , is to note that 

p 7 c 2 + P /max = m / ^( v /n77 + yfnJZ) 2 , (213) 

p/c 2 - P /max = m 7 /i( v /n77 - y/nJI) 2 . (214) 

Thus when P/ max = 0, there are just particles but no antiparticles of the complex 
scalar field, and when P/ max = P/c 2 , there are equal numbers of particles and an- 
tiparticles. This latter possibility is the case when there is only one real scalar field, 
in which case it is a fiction to say that there is the complex scalar field $ at all, but 
the classical real scalar field does act as if it were composed of equal numbers of par- 
ticles and antiparticles of the fictitious complex scalar field, n/+ = nj_ = (l/2)n/. 
(Of course, then there are no vacuum fluctuations of the nonexistent imaginary com- 
ponent of the complex $, but here we are taking the classical limit and are only 
considering the effects of real particles and not of any vacuum fluctuations.) 



By making the arbitrary requirement that the 0(n) transformation lead to real 
ci > C2 > for the coefficients of the two nonzero scalar fields after the transfor- 
mation, we have made an arbitrary choice of what to call particles and what to call 
antiparticles (with the number density of particles never less than the number of 
antiparticles by this choice) and hence of which expression in Eqs. (121 lj) and (|212|) 
has the minus sign in front of the square root. 

Now when we have the possibility of both particles and antiparticles of the 
complex scalar field <&, the U(l) invariance prevents the annihilation of a pair of 
particles or of a pair of antiparticles (at least at the perturbative level) but allows 
the annihilation of a particle-antiparticle pair. On the other hand, the particle- 
antiparticle annihilation cross-section is twice that given by Eq. ()144j) [TH] for two 
real scalar field particles, i.e., 

A-nG 2 m 2 A-nh 2 G 2 m 2 . . 

<r+- = — r- 1 = 7 • 215 

crv c'v 

This means that when we have Nj + particles and Nj_ antiparticles, both of these 
numbers decrease at the rate 

dN I+ diV/_ f A-nh 2 G 2 m 2 1 r rt 
d~t~ = dt~ = J xni + n i- a +- v = 7 J " xn i+ n i- 

tcG 2 r ,„ „ 0> ttG 2 r „ n m 4 , 



c 7 



J d 3 x(2Pf) = ^J d*xP 2 m3X = ^-J d 3 x\^ 2 u \\216) 



The total number decay rate is of course twice this. When m+ = ni-, so that 
-Pfmax = Pic 2 , then the total rate agrees with the integral of Eq. ()147jl . 

When we multiply the total number decay rate for each mass rrij (which is in 
frequency units) by the mass m*j = hmj / c 2 in conventional units and sum over all 
/, we get the total mass loss rate by scalar particle annihilation in conventional mass 
units: 

dM* h 



dt 8vrc 5 



EmS/c^lE^I 2 . (217) 



/ 

We can then multiply this by G/c 3 to get the dimensionless rate of decrease of 
M = GM^/c 3 , the mass in units of time: 

"f = ^X>5/^IX>?,I* (218) 
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13 Quantum and Classical Emission from the 
Simplest Multi-Field Spherical Oscillatons 



Now consider the particular case in which all of the real scalar fields have the same 
mass, so we can drop the subscript I that labels the mass. Furthermore, restrict 
attention to the case in which all of the real scalar fields are oscillating with the 
same quasi-stationary nodeless spherically symmetric mode (that of the simplest 
spherical oscillaton, except that now there are more than one real scalar field that 
may be oscillating with different phases). 

By the argument above, we can perform an 0(n) transformation so that only 
two real scalar fields are then oscillating with nonzero amplitude and are 90° out 
of phase. By the assumption that all of the scalar fields are oscillating in the same 
mode (up to phase), this 0(n) transformation is constant over space, and the ratio 
of the amplitudes of the two resulting nonzero modes are also constant. By the 
procedure above, it can be replaced by a single complex field. 

Let N + be the number of particles of the complex scalar field, N_ be the number 
of antiparticles, and iV = N + + N- be the total number of particles and antiparticles. 
Then if the ip given by Eq. ()76|) (real in this case) represents the simplest spherical 
oscillaton with one real scalar field described above, the complex oscillaton with two 
real scalar fields and the same total mass is represented (after a shift in the origin 
of time) by 

= v (219) 

^2 = -i ^ - V- (220) 

Then the complex scalar field is given by 
1 , 

$ = —=(6i + id> 2 




Now, by using Eq. (j216|) and doing an analysis analogous to that which led to Eq. 
(|152j) . one can show that the quantum annihilation of scalaron particle-antiparticle 
pairs for the spherically-symmetric nodeless complex oscillaton gives 



dt dt 

70 



2Qt^m Li N d N + N„, (222) 



where the numerical constant Q is given in Eq. (jl55j) . 

These equations have the obvious constant of motion being the number of par- 
ticles minus the number of antiparticles, say 

N = N + -N_. (223) 

Then, since the total number of particles is N = N + + iV_, one can write 

N + = ^(N + N ), (224) 

N_ = ±(N-N ). (225) 

Similarly, one can write the dimensionless mass parameter of the classical con- 
figuration as 

/i = Mm = tl x m 2 N (226) 
and also define a constant dimensionless mass parameter as 

/i = t 2 Pl m 2 N . (227) 

As both N + and N_ decay away at equal rates, the oscillaton asymptotically 
approaches the configuration with N = N + = N and AL = 0, which is a static 
boson star with \i = fi . Thus /io is the asymptotic (minimum) value of /i. 

In terms of iV and N , Eq. ()222|) gives 



_ ^ = Qt™m ll N\N 2 - N 2 ) (228) 

from the annihilation of scalaron particle-antiparticle pairs into pairs of gravitons 
(i.e., ignoring the classical emission into scalar radiation). 

Alternatively, we can write the evolution of the dimensionless mass parameter 

as 

-f = Q&™W-ti) ( 229 ) 

from the quantum annihilation. This differential equation has the algebraic solution 

,2 



n 2 



where to is an arbitrary constant of integration. Thus at late times, fi approaches Ho 
exponentially rapidly, and the configuration approaches that of a static boson star 
with the conserved number of particles. 
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It may also be of interest to give the decay rates from the classical emission of 
scalar radiation in this multi-field case with the simplest spherical configuration. By 
using Eq. (JHHj), one can deduce that the classical scalar field emission leads to 

_dN ± = _dN_ _ 2C ^N+N- e „ am 2 i/m 2 N mi) 
dt dt ~ m 5 iV 4 ' K ' 

dN m 6 P1 N 2 -N 2 iJm 2 N 



and 



" I - Cm^e^. (233) 

To get the totals for these rates, one must add the corresponding expressions for 
the quantum rates from Eqs. ()222j) . ()228|) . and ()229j) respectively. For example, the 
total rate at which the mass (in time units) decreases is 

^e-^r + Q&m'M 3 ) (M 2 - Mi), (234) 

where Mq = fio/m = tj^miVo is the asymptotic mass of the final boson star in 
time units, and where the numerical constants C ~ 3 797438, a ~ 39.4338, and 
Q ~ 0.008 513 224 were given to 36 decimal places (though only the first 30-35 of 
these are likely to be accurate) in Eqs. (j!37|) . (|136|) . and (|155J) respectively. 

One can see that at late times, for /i > 0, as /i approaches very near to fio, 
H — fio <C /io/ a > Z 1 approaches /io exponentially rapidly: 

H ~ hq + exp [-(C>o 3 e- a/Ato + Qt 2 Pl m 2 fj$)2mi\. (235) 

Finally, if we use the scalar field mass in conventional mass units, m* = hm/c 2 , 
and the oscillaton mass also in conventional mass units, M* = c 3 M/G, then the to- 
tal mass loss rate from the simplest spherical multiple-field oscillaton with minimum 
conventional mass M*o = c 3 M /G = iV m*, form, < M* < 1O- 1O M [(1 eV)/(m*c 2 )], 
is 

" iJ t » ( c l?^-" , " <Gm • M " , + <«? " <»■ 

Again may I remind the reader that the conventional mass m* of the scalar field is a 
quantum quantity with fi in it, which is why there is an explicit h 2 in the numerator 



72 



of the first (classical) term for the mass decay rate (from scalar field radiation), to 
cancel the implicit Ti 2 in the m 2 term in the denominator. Similarly, in the second 
(quantum) term, once the implicit factor of h 7 in the ml term is taken into account, 
there is one positive power of % appearing, as one would expect for this first-order 
quantum perturbative contribution (the annihilation of pairs of scalarons into pairs 
of gravitons). 

14 Conclusions 

Oscillatons without a U(l) invariance (those without a static geometry) are unstable 
both classically (to emitting scalar waves) and quantum mechanically (to having 
scalarons annihilate into gravitons) . The classical rate dominates for large /x = Mm 
but drops very fast with decreasing /x and is nonanalytic at /x = 0: P c ~ (C//x 2 )e _Q///i , 
Eq. ()135j) . with the numerical constants C given by Eq. ()137j) and a given by 
Eq. (I136j) . The quantum rate also drops as /x drops, but only as a power law in 
/x: P q ~ Q(m/mp\) 2 fi 5 , Eq. ()152j) . with the numerical constant Q given by Eq. 
f!155|) . The quantum rate dominates for /x ~ 1/8 for m ~ H (the current Hubble 
expansion rate, a lower bound on m for any oscillaton existing in our universe today). 

An oscillaton that starts at /Xi « u max w 0.607 jH] has a significant drop in /x 
(more than 10%) over a lifetime comparable with the age of the universe if m*c 2 ~ 
2 x 10~ n eV or M* ~ 3.6 M . However, unless m*c 2 ~ 2.3 x 10 13 eV = 2300 GeV 
or M* ~ 1.8 x 1O~ 24 M « 3.5 x 10 9 g, in this decay time /x does not decrease by 
more than a factor of 2. 

These numerical approximations are using formulas derived in this paper for 
nearly-Newtonian configurations, which have < 1. For accurate results for the 
not-too-small values of /x that would arise from the decay, within astronomical times, 
of oscillatons that started with the maximum mass possible for any reasonable value 
of the scalaron mass, one would need to extend the results derived here to the 
strong-gravity regime. This is research that shall be left to the future. 

Multi-field oscillatons [7j, in which different real scalar fields of the same mass 
oscillate out of phase, do not decay away completely but instead asymptotically 
approach a stable £7(l)-invariant configuration with a static metric (a boson star), 
at a rate given by Eq. jZgp or fl^BEJ). 
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